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ABSTRACT 

Aims. The extensive observations of the supernova SN 1993J at radio wavelengths make this object a unique target for the study of 
particle acceleration in a supernova shock. 

Methods. To describe the radio synchrotron emission we use a model that couples a semianalytic description of nonlinear diffusive 
shock acceleration with self-similar solutions for the hydrodynamics of the supernova expansion. The synchrotron emission, which 
is assumed to be produced by relativistic electrons propagating in the postshock plasma, is worked out from radiative transfer calcu- 
lations that include the process of synchrotron self-absorption. The model is applied to explain the morphology of the radio emission 
deduced from high-resolution VLBI imaging observations and the measured time evolution of the total flux density at six frequencies. 
Results. Both the light curves and the morphology of the radio emission indicate that the magnetic field was strongly amplified in the 
blast wave region shortly after the explosion, possibly via the nonresonant regime of the cosmic-ray streaming instability operating in 
the shock precursor. The amplified magnetic field immediately upstream from the subshock is determined to be B u as 50(//l day) -1 G. 
The turbulent magnetic field was not damped behind the shock but carried along by the plasma flow in the downstream region. Cosmic- 
ray protons were efficiently produced by diffusive shock acceleration at the blast wave. We find that during the first ~8.5 years after 
the explosion, about 19% of the total energy processed by the forward shock was converted to cosmic-ray energy. However, the shock 
remained weakly modified by the cosmic-ray pressure. The high magnetic field amplification implies that protons were rapidly ac- 
celerated to energies well above 10 15 eV. The results obtained for this supernova support the scenario that massive stars exploding 
into their former stellar wind are a major source of Galactic cosmic-rays of energies above ~10 15 eV. We also calculate the flux from 
SN 1993J of gamma-rays arising from collisions of accelerated cosmic rays with ambient material and the result suggests that type 
II supernovae could be detected in 7r°-decay gamma-rays with the Fermi Gamma-ray Space Telescope out to a maximum distance of 
only ~ 1 Mpc. 

Key words, acceleration of particles - magnetic fields - radiation mechanisms: nonthermal - supernovae: individual (SN 1993J) 
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0^ , 1 ■ Introduction is amplified by the diffusive streaming of accelerated particles in 

the upstream region of the shock, which could cause a plasma 

• . Galactic cosmic rays are widely believed to be accelerated in instabilit generating large-amplitude magnetic turbulence (Bell 

. £ ; expanding shock waves initiated by supernova (SN) explosions, & Lucefc 2001 . AmatQ & Blasi 2006; vladimirov et al . 20m 

V at least up to the knee energy of the cosmic-ray spectrum at Magnetic field am pi iflcation via this nonlinear process could fa- 

■ ~10 15 eV and possibly up to the "ankle" at ~10 Iis eV (Axford ciUtate the acceleration of tons in SNRs up t0 „ 10 is eV (e . g . 

CZ | I1994K The theory of diffusive shock acceleration (DSA) of cos- p ar j zot et a \ p2006) 
mic rays is well established (see, e.g., Jones & Ellison [T99U 

Malkov & Drury |200l]for reviews), but two fundamental ques- The acceleration efficiency depends on the diffusive trans- 

tions remain partly unanswered by the theory: what is the max- P ort of energetic particles at both sides of the shock in the self- 

imum kinetic energy achieved by particles accelerated in SN generated turbulence (e.g. Malkov & Drury [200D. The DSA 

shocks and what is the acceleration efficiency, i.e. the fraction theor y cannot accurately predict at the present time how many 

of the total SN energy converted to cosmic-ray energy? particles are injected into the acceleration process as a function 

The maximum cosmic-ray energy attainable in a SN shock of the shock Parameters. Observations of SNRs in high-energy 

scales as the product of the shock size and the strength of the tur- gamma-rays can provide valuable information on the efficiency 

bulent magnetic field in the acceleration region (e.g. Marcowith of cosmic-ray acceleration in these objects. But it remains un- 

et al. 20061. Multiwavelength observations of young shell-type clear whether the emlsslon detected from several SNRs wlth 

supernova remnants (SNRs) provide evidence that the magnetic ground-based atmospheric Cherenkov telescopes is due mainly 

field in the blast wave region is much higher than the interstellar t0 P lon deca y following hadromc collisions of accelerated ions 

magnetic field (e.g. Volk et al.EQ05j Cassam-Chenai et al. 2007 ; or l ° lnverse Compt ° n scatterm S °I^ Wnt phot0nS by aCCeler " 

Uchiyama et al. l2TjQ7l and references therein). A possible expla- ated electrons (e.g. Morlino et al. |2008Jl. 

nation for these observations is that the ambient magnetic field However, theory predicts that efficient acceleration of 

cosmic -ray ions (mainly protons) modifies the shock structure 

Send offprint requests to: V. Tatischeff with respect to the case with no acceleration (see, e.g., Berezhko 

* Permanent address. & Ellison |T9991 >. In particular, the total compression ratio of a 
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cosmic -ray-modified shock can be much higher than that of a 
test-particle shock (i.e. when the accelerated particles have no 
influence on the shock structure), because relativistic particles 
produce less pressure for a given energy density than do nonrel- 
ativistic particles. In addition, the energy loss due to escape of 
accelerated cosmic rays from the shock system can further in- 
crease the compressibility of the shocked gas (see Decourchelle 
et al. 120001) . Furthermore, since the energy going into relativis- 
tic particles is drawn from the shocked-heated thermal popula- 
tion, the postshock temperature of a cosmic-ray-modified shock 
can be much lower than the test-particle value. Observations of 
these nonlinear effects (Hughes et al. 2000; Decourchelle 2005; 
Warren et al. 2005 ) provide indirect evidence for the efficient ac- 
celeration of cosmic -ray protons (which carry most of the total 
nonthermal particle pressure) in shock waves of young SNRs. 

In this paper, we study the production of cosmic-rays by 
nonlinear DSA and the associated magnetic field amplification 
in a very young SN shock, which expands in a relatively dense 
stellar wind lost by the progenitor star prior to explosion. We 
use radio monitoring observations of SN 1993 J conducted with 
the Very Large Array and several other radio telescopes since 
the SN outburst (see Weiler et al. [20071 Bartel et al. I20071 and 
references therein). The radio emission from SNe is thought to 
be synchrotron radiation from relativistic electrons of energies 
< 1 GeV accelerated at the expanding blast wave (Chevalier 
1982b; Fransson & Bi6rnsson |19981 ). The key motivation for the 
present study is that the radiating electrons can be influenced by 
the presence of otherwise unseen shock-accelerated protons. In 
particular, the theory of nonlinear DSA predicts that the energy 
distribution of nonthermal electrons below 1 GeV steepens with 
increasing efficiency of proton acceleration (Ellison et al. 2000). 

We chose SN 1993J because it is one of the brightest radio 
SNe ever detected and has already been the subject of numerous 
very useful observational and theoretical studies (see, e.g., for 
the radio emis sion of SN 1993 J B artel et al. rT994l [20001 [20071 
Marc aide et al. [19941 [T9951 [19971 Van Dyk et al. [19941 Fransson 
et al. 119961 Fransson & Bjornsson 1998; Perez-Torres et al. 
I200T1 Miod uszewski et al. [20011 Bietenholz et al. [20031 Weiler 
et al.[2007 and references therein). 

The model of the present paper is largely based on the work 
of Cassam-Chenai et al. (2005; see also Decourchelle et al. 2000; 
Ellison & Cassam-Chenai 2005 ) on the morphology of syn- 
chrotron emission in Galactic SNRs. These authors have devel- 
oped a model to study radio and X-ray images of SNRs undergo- 
ing efficient cosmic -ray production, that couples a semianalytic 
description of nonlinear DSA with self-similar solutions for the 
hydrodynamics of the SN expansion. The model has recently 
been applied to the remnants of Tycho's SN (Cassam-Chenai' et 
al. 120071) and SN 1006 (Cassam-Chenai et al. l2008l . 

Duffy et al. (119951 ) have studied the radio emission from 
SN 1987 A by considering a two-fluid system consisting of a 
cosmic -ray gas and a thermal plasma to calculate the structure of 
the blast wave. In their model, however, the efficiency of cosmic- 
ray acceleration is not deduced from the radio data, but assumed 
to be similar to that required in SNRs to explain the observed 
flux of Galactic cosmic rays. Another difference with the present 
work is that the turbulent magnetic field in the shock precursor 
is not assumed to be amplified by the DSA process, but is taken 
by these authors to be of the same order as the ordered field in 
the wind of the progenitor star. 

A preliminary account of the present work has been given 
elsewhere (Tatischeff 2008 ) and all of the present results super- 
sede those published earlier. 



2. Model 

2.1. Forward shock expansion 

The type lib SN 1993 J was discovered in the galaxy M81 by 
Garcia (Ripero et al. [19931 o n 199 3 March 28, shortly after 
shock breakout (Wheeler et al. 119931 ). Very long baseline inter- 
ferometry (VLB I) observations (see, e.g., Marcaide et al. 1997; 
Bietenholz et al. 12003 ) revealed a decelerating expansion of a 
shell-like radio source. The radio emission is presumably pro- 
duced between the forward shock propagating into the circum- 
stellar medium (CSM) and the reverse shock running into the SN 
ejecta (e.g. Bartel et al. 120071) . Marcaide et al. ( 119971) found the 
outer angular radius of the radio shell to evolve as a function of 
time t after explosion with the power law O cc t m , where the de- 
celeration parameter m = 0.86 + 0.02. More recently, Weiler et 
al. ([2007 and references therein) reported that the angular expan- 
sion of SN 1993 J up to day 1500 after outburst can be expressed 
as 6 = 6.2 x (t/l day)'" //as with m = 0.845 ± 0.005; here is by 
definition the angular radius of the circle that encompasses half 
of the total radio flux density to better than 20%. 

These results are consistent with the standard, analytical 
model for the expansion of a SN into a CSM (Chevalier 1982a 
1983; Nadyozhin 1985). This model assumes power-law den- 
sity profiles for both the outer SN ejecta, p e j — C^f^R^" (with 
n > 5), and the CSM, p C SM = CiR s (with s < 3), where C 2 
and C\ are constants. The SN expansion is then found to be self- 
similar (i.e. the structure of the interaction region between the 
forward and reverse shocks remains constant in time except for a 
scaling factor) and the deceleration parameter m - (n-3)/(n-s). 
For a standard wind density profile with s = 2 (see below), 
the deceleration reported by Weiler et al. (J2007 ) corresponds to 
n w 8.5, in fair agreement with numerical computations of SN 
explosions (e.g. Arnett [T988l . 

However, Bartel et al. (2000 2002) found significant 
changes with time of the parameter m, indicating deviations from 
a self-similar expansion. These authors determined the outer an- 
gular radius o as a function of time by consistently fitting to the 
two-dimensional radio images observed at 34 epochs between 
1993 and 2001 the projection of a three-dimensional spherical 
shell of uniform volume emissivity. They fixed the ratio of the 
outer to inner angular radius (as expected from a self-similar ex- 
pansion) at 6 /6i = 1.25. By performing a least-squares power- 
law fit to the values of 6 thus determined (see Fig.[T|), they ob- 
tained a minimum reduced x 1 of xl - 1 -8 for 64 degrees of 
freedom, indicating that the null hypothesis of a self-similar ex- 
pansion can be rejected. 

As discussed by Bartel et al. (2002 ), this result strongly de- 
pends on the systematic errors associated with the model of uni- 
form emissivity in the spherical shell. The overall systematic un- 
certainty was estimated to range from 3% at early epochs to 1% 
at later epochs and to dominate the statistical error in most of the 
cases. We show below (Sect. 3.1) that the assumption of uniform 
emissivity is questionable, because radio-emitting electrons ac- 
celerated at the forward shock are expected to lose most of their 
kinetic energy by radiative losses before reaching the contact 
discontinuity between shocked ejecta and shocked CSM (see 
also Fransson & Bjornsson 1998). Furthermore, the radial emis- 
sivity profile is expected to vary with time and radio frequencjQ. 
Although the systematic uncertainties were carefully studied by 
Bartel et al. (2002), we note that a moderate increase of the er- 



1 This frequency dependence may explain why the outer radii ob- 
tained by Bartel et al. (2002 ) at 2.3 and 1.7 GHz tend to be larger by up 
to 5% than those measured by the same authors at higher frequencies. 
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Fig. 1. Time evolution of the outer angular radius of the shell- 
like radio emission from SN 1993J. The data were obtained by 
Bartel et al. ( [20021 ) from observations at 1.7, 2.3, 5.0, 8.4, 14.8, 
and 22.2 GHz. The dotted line is a least-squares power-law fit to 
the data (Eq.Q}. 



rors would make the expansion compatible with the self-similar 
assumption. For example, by setting a lower limit of 3% on the 
uncertainties in the outer radii measured by Bartel et al., we get 
from a least-squares power-law fit to the data xl — 1 -25 and an 
associated probability of chance coincidence of 8%. Thus, the 
null hypothesis of a self-similar expansion could not anymore 
be rejected at the usual significance level of 5%. 

For simplicity, I shall use the self-similar solution to model 
the hydrodynamic evolution of the SNR (Sect. 2.5). The radius 
of the forward shock, R s — 6 D where the source distance D = 
3.63 + 0.34 Mpc (Freedman et al. |1994| l, is estimated from a 
power-law fit to the data of Bartel et al. (2002 ) at all frequencies 
(Fig-EK which gives 



G = (0.292+0.004) 



100 days 



0.829±0.005 



mas. 



(1) 



Thus, we have R s = R (t/1 day)'" with R = 3.49 x 10 14 cm and 
m = 0.83. The forward shock velocity is then 



v _ dRs _ y I f 
s dt \ 1 day 



with Vq = 3.35 x 10 4 km s" 1 . (2) 



In the self-similar model, the forward shock radius is given 
by (e.g. Chevalier [T983l 



'Ci 



!/(«-*) 



(3) 



where the constant K s depends on s, n, and the adiabatic index 
of the shocked gas y g . For s — 2, n — 7.88, corresponding to 
m = (n - 3)/(n - s) = 0.83, and y g = 5/3 (see below), we 
have K s = 0.93. The parameter C\ that fixes the CSM density 
can be estimated from the radio emission. Anticipating the re- 
sults presented in Sect. 3, the best-fit model to the radio light 
curves gives C\ = 1.9 x 10 14 g cm -1 . Comparison of Eq. (01 
with the observed SN expansion then provides C2, which in turn 
can be expressed in terms of the explosion energy and ejected 
mass (e.g. Nadvozhin |1985t Decourchelle & Ballet |1994l l. We 
finally obtain 



where / e j is a parameter of order unity that depends on the veloc- 
ity distribution of the inner ejecta and £51 is the kinetic energy 
of the explosion in units of 10 51 erg. This result is consistent 
with the ejected mass estimated from the optical light curve of 
SN 1993J, 1.9 < M ej < 3.5 M (Young et al. [19951 . However, 
numerical simulations of SN explosions show that the density 
profile of the outer ejecta can be more complicated than a power 
law (see, e.g., Iwamoto et al. 1997) and one should bear in mind 
that the self-similar solution is a strong simplification of the real 
hydrodynamic evolution. 



2.2. Free-free absorption of the radio emission in the 
circumstellar medium 

Figure [2] shows a set of light curves measured for SN 1993 J at 
0.3 cm (85-110 GHz), 1.2 cm (22.5 GHz), 2 cm (14.9 GHz), 
3.6 cm (8.4 GHz), 6 cm (4.9 GHz), and 20 cm (1.4 GHz). We see 
that at each wavelength the flux density first rapidly increases 
and then declines more slowly as a power in time (the data at 
0.3 cm do not allow to clearly identify this behavior). The ra- 
dio emission was observed to suddenly decline after day ~3100 
(not shown in Fig. [2J, which is interpreted in terms of an abrupt 
decrease of the CSM density at radial distance from the progen- 
itor R out ~ 3 x 10 17 cm (Weiler et al. 120071) . This is presumably 
the outer limit of the dense cocoon which was established by a 
high mass loss from the red supergiant progenitor of the SN for 
~10 4 years before explosion. 

We see in Figure [2] that the maximum intensity is reached 
first at lower wavelengths and later at higher wavelengths, which 
is characteristic of absorption processes. For SN 1993J, both 
free-free absorption (FFA) in the CSM and synchrotron self- 
absorption (SSA) are important (Chevalier 119981 Fransson & 
Bjornsson [19981 Weiler et al. 120071 . Th e Razin effect can be 
excluded (Fransson & Bjornsson 119981 ). FFA is produced in 
the stellar wind that has been heated and ionized by radiation 
from the shock breakout. If the circumstellar gas is homoge- 
neous and of uniform temperature 7csm, the radio emission pro- 
duced behind the forward shock is typically attenuated by a fac- 



tor exp(- 



homog 
r CSM 



satisfies (Weiler et al. 1 19861 ) 



), where the optical depth at a given frequency v 



homog 
r CSM 



Jr*oo 
2 T -1.35 -2.1 , p „ tt-1,35, -2.1pl-2s 
Pcsm^csm 17 dRccl n ^„V R, 



CSM 



oc T, 



•1.35-2.1 
CSM y 



jn(l-2j) 



(5) 



(4) 



From a fit to early data using the parametrized model of Weiler 
et al. (119861120021 see Appendix A), Van Dyk et al. JT9941 ) found 
for the time dependance of the optical depth 6 = m(l - 2s) ~ -2. 
They concluded that the density profile of the CSM must be sig- 
nificantly flatter, s ~ 1.5, than that, s — 2, produced by a con- 
stant mass-loss rate and constant-velocity stellar wind. They in- 
terpreted this result in terms of a steady decrease of the mass-loss 
rate of the SN progenitor star prior to explosion. This conclusion 
was later confirmed by Fransson et al. (1996). But in both studies 
SSA was not taken into account. 

Fransson & Bjornsson (1998 ) performed the most detailed 
modeling of the radio emission from SN 1993J to date. In par- 
ticular, they took into account both FFA and SSA, and included 
all relevant energy loss mechanisms for the relativistic electrons. 
As a result, they were able to adequately reproduce the radio 
light curves with the standard s - 2 density profile. In their 
model, the measured time dependence of the FFA optical depth 
(i.e. 5) is accounted for by a decrease of Tcsm with radius like 
T^csm K R~ l (see also Fransson et al. 1 19961 ). However, as pointed 



4 



V. TatischefF: Radio emission and nonlinear diffusive shock acceleration in SN 1993J 



10 
10 



10 



a 

3 



10 



10 



. 1 1 1 1 ' ' '. 

0.3 cm _ 


■ i i i ■ 1 ■ '. 

1.2 cm 

r / / ^^*h£± i 


\ i 1 1 1 ni| — i i 1 1 1 ni| — i i 1 1 ini| — i 
2 cm 


\ 1 1 1 mil — i 1 1 1 1 mi — i 1 1 1 1 mi — i 1 1 \. 

'. ••jafifc. ^ cm - 




''III llll| 1 1 1 1 1 lll| 1 1 1 1 1 lll| Mil 

■ 20 cm 



10 10 10 

Days after outburst 



10 10 10 
Days after outburst 



Fig. 2. Radio light curves for SN 1993J at 0.3, 1 .2, 2, 3.6, 6, and 20 cm. The data are from Weiler et al. (2007) and references therein. 
The solid lines represent the best-fit semi-empirical model of these authors. The dotted lines show a synchrotron self-absorption 
(SSA) model in which the SSA optical depth and the unabsorbed flux density were obtained from fits to radio spectra taken between 
day 75 and day 923 after outburst (see Sect. 2.3). The difference between the latter model and the data before day 100 is partly due 
to free-free absorption of the radio emission in the CSM. 



out by Weiler et al. (2002), no evidence for such a radial depen- 
dence of Tqsm is found in other radio SNe (e.g. SN 1979C and 
SN 1980K). 

Immler et al. ((2001) provided support from X-ray observa- 
tions to the scenario of a flatter CSM density profile, as they 
found s = 1.63 from modeling of the observed X-ray light curve. 
Their analysis assumes that the X-ray emission arises from the 
forward, circumstellar shock. But later than ~200 days post- 
outburst, the X-ray radiation was more likely produced in the 
SN ejecta heated by the reverse shock (Fransson & Bjornsson 
12005 1 and references therein). 

With the model of the present paper, the measured radio light 
curves can be well explained with the standard s = 2 assumption, 
but not with a much flatter CSM density profile. This result is 
independent of the admittedly uncertain FFA modeling, as we 
will see in Sect. 4. 1 that with s — 1.6 the optically thin emission 
cannot be simultaneously reproduced at all wavelengths in the 
framework of the model (see Fig. [14]). So the question is how a 
Pcsm K R 2 density profile can be reconciled with the relatively 
low value of 6 implied by the data. 

A possible explanation is that the absorbing CSM is inhomo- 
geneous. Both radio (Weiler et al. 1990) and optical and UV (e.g. 
Tran et al. 119971 ) observations show that the wind material lost 
from SN progenitors can be clumpy and/or filamentary. FFA of 
the radio emission by a nonuniform CSM can be accounted for 
by an attenuation factor of the form [1 - e x P( - '*"pg^ ps )]/Tp5^ p \ 

where T^™ ps is the maximum of the optical depth distribution, CSM _ \Anf c \u w nin 



which depends on the number density and geometric cross sec- 
tion of the clumps of wind material (Natta & Panagia [1984; 
Weiler et al. 2002). Weiler et al. (2007 ) performed an overall 
fit to all of the measured radio light curves of SN 1993 J using a 
parametrized model that takes into account both SSA and FFA, 
and includes both attenuation by a homogeneous and inhomo- 
geneous CSM (see Appendix A). In their best-fit model (solid 
lines in Fig. [2]), FFA is mainly due to the clumpy CSM, the cor- 
responding optical depth, T^™ ps , being much larger than t^™ 8 
in the optically thick phase for all wavelengths (see Table 4 in 
Weiler et al. 11507]). 

In view of these results, we are going to neglect in first ap- 
proximation the attenuation of the radio emission by the homo- 
geneous component of the CSM in front of the attenuation by 
the clumpy CSM, which will allow us to derive in a simple way 
and self-consistently the mass loss rate of the progenitor star and 
the structure of the radio emission region (see below). We antici- 
pate, however, that the adopted simple FFA model will not allow 
us to accurately reproduce the rising branches of the light curves 
in the optically thick phase. 

Following Weiler et al. ( 1986 2002), the optical depth pro- 
duced by FFA in a clumpy presupernova wind can be written 

as 



clumps 



RSG 



1 +2X\(k { - { 



1+4XJ\3R 3 Q 



1 day 



(6) 
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with, for a uniform 7csm, 



Kf _ { = 6.34 x 10 



-29 



v-2.1/ rj, \ -1.35 

cm' . (7) 



5 GHz/ \2xl0 5 K 



Here, Mrsg and u w are, respectively, the mass loss rate and wind 
terminal velocity of the red supergiant progenitor, ntjj is the mass 
of a hydrogen atom, X is the He to H abundance ratio in the pre- 
supernova wind, and f c \ is a factor that depends on the number 
and geometrical properties of the clumps. Weiler et al. (2002) 
argue that f c \ * 0.67 (resp. / c i ~ 0.16) for attenuation by a sta- 
tistically large (resp. small) number of clumps along the line of 
sight. By equating t c ( ^^ s from Eq (0 with its expression in the 
parametric model of Weiler et al. summarized in Appendix A 
(Eq. IA.61) . one can derive a relation for the progenitor mass loss 
rate as a function of the normalization parameter ^3: 



M RSG = 2.2 x 10" 7 V#^)(: 



?CSM 

lOkms'/UxlOSK 



0.675 



Ro 



,3.49 x 10 14 cm 



1.5 



M yr 



(8) 



where we adopted X = 0.3 (e.g. Shigeyama et al. |19941 >. In the 
following we use u w = 10 km s _1 and Tqsm = 2xl0 5 K, the latter 
value being based on the photoionization calculations performed 
for SN 1993J by Fransson et al. (I1996K We note, however, that 
according to these calculations 7csm was higher during the first 
~10-20 days post-outburst (see Fig. 1 1 of Fransson et al. 1996). 

FFA in the progenitor wind of SN 1993 J is estimated below 
from the following iterative process. First the synchrotron emis- 
sion as a function of time after outburst is calculated from a set 
of reasonable initial parameters. External FFA is then estimated 
from an overall fit of the calculated light curves to the radio data, 
using Eqs. (1A.3I ) and dA.6t for the CSM attenuation factor with 
fixed 5' - -3m (see Eq.|6]l. Thus, the only free parameter is ^3. 
The corresponding progenitor mass loss rate is obtained from 
Eq. ((S). The derived wind density just upstream from the for- 
ward shock, p„ = Mrsg/(4tt/?jM vv ), is then used to calculate the 
shock properties and the associated synchrotron emission. The 
process is continued until convergence is reached. 



2.3. Time evolution of the mean magnetic field in the 
synchrotron-emitting region 



In the original model of Chevalier (1982b) for the radio emis- 
sion from SNe, the magnetic energy density in the radio-emitting 
shell is assumed to scale as the total postshock energy density 
(oc p„V 2 ), such that the strength of the mean magnetic field 
(B) oc r 1 for s = 2. Later on, Chevalier (fl996l [T9981 also con- 
sidered that the postshock magnetic field could result from the 
compression of the circumstellar magnetic field, which would 
imply (B) oc R~ l oc r m . Because SSA plays an important role in 
the radio emission from SN 1993 J, the evolution of the mag- 
netic field can be estimated from the measured light curves. 
Fransson & Bjornsson (1998) found, however, the two scaling 
laws (B) oc f _I and (B) oc R~ l to be compatible with the data for 
SN 1993 J available at that time. 

We re-estimate here the temporal evolution of {B) using the 
fitting model of Weiler et al. ( 19861 120021 ) together with the sim- 
ple formalism given in Appendix A. We analyze individual radio 
spectra at various dates of observation to determine the time vari- 
ation of the fitted parameters. The method is similar to the one 
previously employed by Fransson & Bjornsson ( 119981 ). 



We use observations made with the Very Large Array (VLA) 
in which a non-zero flux density was measured at at least five 
wavelengths at the same time (1.2, 2, 3.6, 6, and 20 cm). The cor- 
responding data were obtained later than 75 days after explosion 
(see Fig. [2](. We also use the data at 90 cm taken at day 922.7 
post-outburst. Each radio spectrum is fitted with the function 

F v = Ki(v/5 GHz)°'Apg™ ps 'AssA, where the attenuation factors 
are given by Eqs. dA.3t and iAAi , with the corresponding optical 
depths Tcsm PS = £s(v/5 GHz)" 21 and t SS a = K 5 (v/5 GHz)" 25 . 
The fitting function thus contains four free parameters: K\, a, 
£3, and £5. 

The best-fit parameters are shown in Figure [3] We took into 
account only the spectral fits of relatively good quality, with^ 2 < 
2. It was checked that the final result is not strongly dependent 
on this data selection. The FFA optical depth was found to 
be compatible with zero at most epochs. No regular evolution of 
this parameter can be deduced from the fitting results (Fig. [3};). 
On the other hand, the temporal evolutions of K\ and K5 can be 
well described by power-law fits (Figs. [3^ and[3ji). We found 



Ki = (185.3 + 5.1) 



100 days 



K 5 = (1.253 + 0.090) 



100 days 



mJy, with B = -0.731 + 0.026 



, with 5" = -2.971 + 0.094, 



(9) 



with reduced^- 2 of 1.2 and 1.3, respectively. The spectral index 
a does not show a significant evolution with time (Fig.rjji). Its 
average constant value is a — -0.887 + 0.020 (xl = 0.32). Then, 
from Eq. iAA% with y = 1 - 2a = 2.11 A + 0.040, we get 



(B) = (2.4+1.0) 



100 days 



G, with b = -1.16 + 0.20. (10) 



We have neglected here the correlations between the various fit- 
ted parameters in the error determination, which is a good ap- 
proximation given that the error in the magnetic field at day 100 
mainly arises from the uncertainty in the source distance (9.4%) 
and the error in the power-law index b is dominated by the un- 
certainty in 6" . 

The mean magnetic field given by Eq. ( TTOb is in good agree- 
ment with the one obtained by Fransson & Bjornsson (1998). 
But in contrast with the conclusions of these authors, the value 
of b obtained in the present analysis indicates that the scaling law 
{B) oc R' 1 can be excluded at the 90% confidence level. The dis- 
crepancy is partly due to the different assumptions made for the 
SN expansion. Indeed, Fransson & Bjornsson (1998) assumed 
m = 1 for t < 100 days (i.e. R s oc f) and m = 0.74 at later epochs. 
The FFA model is also different in the two analyses. However, 
because the present result is restricted to relatively late epochs, 
it is only weakly dependent of the FFA modeling. 

A "pure" SSA model calculated with a = -0.887 and the 
best-fit power laws for ^1 and ^5 (Eq.|9]( is shown in Figure [2] 
We see that this model does not correctly represent the data taken 
before day 100. This is partly because free-free attenuation of 
the radio emission in the CSM was not taken into account. But 
it is also due to the strong radiative losses suffered by the radio- 
emitting electrons at early epochs (see Fransson & Bjornsson 
1998 and Sect. 3 below), which are not included in the paramet- 
ric formalism of Weiler et al. (2002 and references therein) used 



here. As we will see below, one of the main effects of the electron 
energy losses is to reduce the flux density at short wavelengths 
during the transition from the optically thick to the optically thin 
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Fig. 3. Evolution of the parameters K\ , a, and ^5 obtained from fits to individual VLA spectra (see text). Formally, K\ represents 
the unabsorbed flux density at v — 5 GHz, a the corresponding synchrotron spectral index, K3 and K$ the FFA and SSA optical 
depths, respectively, also at v — 5 GHz. The dotted lines show least-squares fit to the extracted parameters (see text). 



regime. If this effect is not properly taken into account, the rela- 
tive contributions of FFA and SSA at early epochs post-outburst 
cannot be reliably estimated. Consequently, the magnetic field 
determination for these epochs is uncertain. We assume here that 
the time evolution of {B) determined from day 75 after explosion 
was the same at earlier times. 

We also see in Figure [2] that the decline with time of the 
optically-thin emission calculated in the SSA model is too slow 
as compared to the data at 20 cm. A possible explanation is that 
the energy distribution of the radiating electrons deviates from 
a power law of constant spectral index, as implicitly assumed in 
the above modeling (but see below). 

The derived time dependence of (B) is consistent with the 
scaling law originally adopted by Chevalier ( |1982bl l, (B) oc r 1 . 
In Chevalier's model, the scaling law is based on the assump- 
tion that the postshock magnetic field is built up by a turbulent 
amplification powered by the total available postshock energy 
density. We will see in Sect. 4.3 that the temporal evolution 
(B) oc t is also to be expected if the magnetic field is am- 
plified by Bell (2004)'s nonresonant cosmic-ray streaming in- 
stability in the shock precursor region and if furthermore the 
shock is not strongly modified by the back pressure from the 
energetic ions. We note that the time dependence (B) oc r l 
was also recently reported by Soderberg et al. (2008) for the 
Type Ibc SN 2008D. Therefore, based on both observational evi- 
dences and a theoretical basis, I assume in the following that the 
postshock magnetic field results from an amplification by the 
cosmic -ray streaming instability operating in the shock precur- 
sor and that the field immediately upstream from the subshock is 
of the form B„ = B„o(t/l day) -1 , where B u q is a free parameter 
to be determined from fits to the radio data. 

2.4. Nonlinear particle acceleration 

Particle acceleration at the forward shock is calculated with the 
semianalytic model of nonlinear DSA developed by Berezhko 



& Ellison (TT9991 and Ellison et al. d2000l Although the model 
strictly applies to plane-parallel, steady state shocks, it has been 
successfully used by Ellison et al. (200Q for evolving SNRs and 
more recently by Tatischeff & Hernanz (120071) to describe the 
evolution of the blast wave generated in the 2006 outburst of 
the recurrent nova RS Ophiuchi. The main feature of this rel- 
atively simple model is to approximate the nonthermal part of 
the shocked proton and electron phase-space distributions as a 
three-component power law with an exponential cutoff at high 
momenta, 

f P (p) = a p p-^ ) CKp(-p/p p max ) (11) 
and 

f e (p) = a e p-^ exp(-W/Cx) ■ (12) 

where the power-law index q(p), which is the same for pro- 
tons and electrons, can have three different decreasing values 
in the momentum ranges p < m p c (m p is the proton mass and 
c the speed of light), m p c < p < 0.01/?m ax , and p > 0.01p^ ax . 
This piecewise power-law model is intended to account for the 
upward spectral curvature that results from efficient ion accel- 
eration. The number density of particles per unit energy inter- 
val, N(E), is related to the phase-space distribution function by 
N(E) = 4np 2 f(p)dp/dE. 

The maximum proton momentum p p mm is calculated either 
by time integration of the DSA rate (i.e., shock age limitation) 
or by equalling the upstream diffusion length to some frac- 
tion /esc of the shock radius (i.e., particle escape limitation), 
whichever produces the lowest value of /?jj, ax (see, e.g., Baring 
et al. I1999I I. Following, e.g., Ellison & Cassam-Chenai (2005), 
I take /esc = 0.05. To estimate the spatial diffusion coefficient, 
k=Av/3, the scattering mean free path A of all particles of speed 
v is assumed to be A — r) m f p r g (Ellison et al. I2000t . where r g is 
the particle gyroradius and 77 m f p is a constant that characterizes 
the scattering strength. I use ?7m fp = 3 , which is a typical value 
for young SNRs (Parizot et al. 2006). The maximum electron 
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Fig. 4. Shocked proton and electron phase-space distributions vs. kinetic energy, at day 1000 after shock breakout. Following 
Berezhko & Ellison ( 119991 1, the phase-space distribution functions have been multiplied by [p/(m p c)] 4 to flatten the spectra, and 
by [{m p c) 3 /n u ] to make them dimensionless (n u is the proton number density ahead of the shock precursor). The upper curves are 
for protons and the lower ones for electrons. The three sets of injection parameters i/f . and rjf. are those used for the synchrotron 
calculations shown in Fig. Qj] The magnetic field used for these calculations is B u = B„o(f/l day) -1 = 50 mG (i.e. B M = 50 G). 



momentum p e m3X is limited by synchrotron and inverse Compton 
losses (see Sect. 2.6). 

Given the upstream sonic and Alfven Mach numbers of the 
shock, which can be readily calculated from V s , Tqsm, Pu, and 
B u q (see below), the proton distribution function (i.e. the nor- 
malization a p and power-law index q{p)) is determined by an 
arbitrary injection parameter 77?., which is the fraction of to- 
tal shocked protons in protons with momentum p>p'^ injected 
from the postshock thermal pool into the DS A process. The work 
of Blasi et al. (2005) allows us to accurately relate the proton in- 
jection momentum pf . to nf .. 
j * ln j / ln j 

The normalization of the electron distribution function is ob- 
tained from (Ellison et al. |20001 > 



''inj / m e 

a p— 



'inj 



m 



(? S ub-3)/2 



(13) 



where rf^. is the electron injection parameter (i.e. the fraction of 



shocked electrons with superthermal energies), m e is the elec- 
tron mass, and q su h the power-law index q(p) for p < m p c. 
Nonthermal particles in this momentum range are accelerated 
at the gas subshock and we have (Berezhko & Ellison T9991 > 



3r sub 
r su b - 1 



(14) 



where r su b is the compression ratio of the subshock. 

Alfven wave heating of the shock precursor is taken into ac- 
count from the simple formalism given in Berezhko & Ellison 
( 119991 1. However, a small change to the model of these authors 
is adopted here: the Alfven waves are assumed to propagate 
isotropically in the precursor region and not only in the direction 
opposite to the plasma flow, i.e. Eqs. (52) and (53) of Berezhko 
& Ellison ( 1999) are not used. This is a reasonable assumption 
(see, e.g., Bell & Lucek 2001) given the strong, nonlinear mag- 
netic field amplification required to explain the radio emission 
from SN 1993J. 

Figure [4] shows calculated shocked proton and electron 
phase-space distributions for three sets of injection parameters 
(rf ..if .) that will be used in Sect. 3.2 to model the radio light 
curves. The thermal Maxwell-Boltzmann components were cal- 
culated using the shocked proton temperature T P determined 



by the nonlinear DSA model (see, e.g., Ellison et al. 2000) 
and arbitrarily assuming the temperature ratio T e JT p = 0.25. 
Noteworthy, the nonthermal electron distribution is independent 
of T e s /T P when 77^. is specified (see Eq.fTSl except for the elec- 
tron injection momentum 



Pim = P 



m e T e s 
m„T p 



1/2 



(15) 



The uncertain temperature ratio has thus practically no influence 
on the modeled radio emission. 

We see in Figure that for 77^ = 10~ 5 the well-known 

test-particle result q(p) = 4 is recovered. But for 77^. > 10~ 4 
(Figs. [4J3 and c) the nonlinear shock modification becomes sig- 
nificant. In particular we see that the nonthermal electron distri- 
bution steepens below ~ 1 GeV with increasing 77?., as a result 
of the decrease of r su b (Eq.[T4li. This is important because the ra- 
dio emission from SN 1993 J is produced by relativistic electrons 
of energies < 1 GeV. 

Figure [5] shows calculated subshock and total compression 
ratios for the case rf . = 10~ 4 which, as will be shown in 

'inj 

Sect. 3.2, provides the best description of the radio light curves. 
The calculations were performed with the upstream sonic Mach 
number M$ >u = V s /cs )U = 560(f/l day)" 1-1 , given the upstream 
sound velocity cs, u - [Tsj^csm/^'Wh)] 1 ^ 2 = 60 km s _1 for 
Tcsm — 2 x 10 5 K. Here y g = 5/3 is the adiabatic index 
for an ideal non-relativistic gas, k is the Boltzmann's constant, 
and p. = (1 + 4X)/(2 + 3X). Given the magnetic field imme- 
diately upstream from the subshock B u = B„o(f/l day) -1 , the 
Alfven Mach number Ma,u - V s /ca, u is independent of time; 
here ca,u = Bu/ y/4-npu is the Alfven velocity. Anticipating the 
results presented in Sect. 3, with the best-fit parameter values 
B u0 = 50GandM RSG = 3.8xl0~ 5 M yr" 1 , we have M A; „ = 9.5. 
Thus, Ma, u <S u , which implies that energy should be very ef- 
ficiently transferee! from the accelerated particles to the thermal 
gas via Alfven wave dissipation in the shock precursor region 
(Berezhko & Ellison 119991) . The resulting increase in the gas 
pressure ahead of the viscous subshock limits the overall com- 
pression ratio, r tot , to values close to 4 (i.e. the standard value for 
a test-particle strong shock). 

However, we see in Figure [5] that the acceleration efficiency 
e nt increases with time. This quantity is defined as the fraction of 
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as a function of time after outburst, for 77^ = 10~ 4 and the up- 
stream sonic and Alfven Mach numbers of the shock Ms ,„ = 
560(?/l day)"' -1 and M Aj „ = 9.5, respectively (see text). 



total incoming energy flux, Fo = Q.5p u V^, going into shock- 
accelerated nonthermal particles. At day 3100 after outburst, 
when the shock has reached the outer boundary of the dense 
progenitor wind, e nt = 23%. The subshock compression ratio 
is found to slowly decrease with time, as the shock becomes in- 
creasingly modified. Thus, we expect a gradual steepening with 
time of the electron distribution between pf^ and m p c (Eq.[T4"li. 

Given the high turbulent magnetic field in the shock region, 
the timescale for diffusive acceleration of electrons in this mo- 
mentum range is very rapid, r acc «: 1 hour. This is much shorter 
than the characteristic timescale for variation of the shock struc- 
ture (see Fig. 0, so it is justified to assume that the spectrum of 
the radio-emitting electrons at a given time is determined by the 
instantaneous subshock compression ratio at that time. We note, 
however, that our calculation of the high-energy end of the pro- 
ton spectrum is not accurate, because the acceleration timescale 
for the highest energy particles is much longer. 
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Fig. 6. Radial profiles of (a) the gas density and (b) the total 
magnetic field in the shock region. The density is normalized 
to the upstream value p„ = Mrsg/(4-7tR 2 s ii w ). The magnetic field 
is normalized to the immediate postshock value B4 — B u [(l + 
2r t 2 0t )/3]^ 2 . In panel (b), the solid line shows the profile of the 
postshock magnetic field carried by the flow and the dashed lines 
that of the damped magnetic field at days 100 and 1000 after 
outburst (Eqs.|2T|and|27]with Ma,u = 9.5). The magnetic field is 
not modeled inside the contact discontinuity. 



2.5. Magnetohydrodynamic evolution of the postshock 
plasma 

The hydrodynamic evolution of the postshock plasma is cal- 
culated using the self-similar model (Chevalier |1982al [1983; 
Nadyozhin 119851 with the deceleration parameter m = 0.83 
(Sect. 2.1), the standard s — 2 density profile (Sect. 2.2), and 
the adiabatic index y g = 5/3. Thus, the effects of the back 
pressure from the accelerated ions on the dynamics of the SNR 
are neglected. It is a good approximation for SN 1993J given 
that 4 < r tot < 4.4 for the best parameter value 77^ = 10~ 4 
(Fig-El. The situation is different in older SNRs, such as the rem- 
nant of Kepler's (Decourchelle et al.[2000) and Tycho's (Warren 
et al. 120051 1 SNe. In these objects, the backreaction of shock- 
accelerated cosmic rays has more influence on the shock struc- 
ture, mainly because the magnetic field in the precursor region 
is much lower than for SN 1993J, such that Alfven wave heating 
is less important. 

Figure [6^ shows the density profile in the region of inter- 
action between the SN ejecta and the progenitor wind. The lo- 
cations of the forward and reverse shocks are clearly visible, at 



1.27 and 0.976 times the radius of the contact discontinuity, re- 
spectively. 

The postshock magnetic field is thought to result from the 
shock compression of the turbulent magnetic field immediately 
upstream from the subshock, that has been presumably ampli- 
fied in the precursor region by both resonant (Bell & Lucek 
2001) and nonresonant (Bell 2004) cosmic-ray streaming insta- 
bilities. Following the work of Cassam-Chenai' et al. (20Q7J for 
Tycho's SNR, I make two different assumptions about the post- 
shock magnetic field evolution: one where the turbulent mag- 
netic field is simply carried by the downstream plasma flow (i.e. 
advected) and another where the magnetic turbulence is rapidly 
damped behind the blast wave. 

2.5.1. Advected downstream magnetic field 

The equations describing the evolution of the postshock mag- 
netic field in the case of pure advection in the downstream 
plasma are given in Cassam-Chenai et al. ((2005) and references 
therein. They are reproduced here for sake of convenience. Let 
the radial and tangential components of the immediate postshock 
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magnetic field be Bd. r and Bdj, respectively. Assuming that the 
upstream magnetic field is fully turbulent and isotropic, we have 



B A , = 1/V3B„ 
B d j = ^J2/3r tot B u 



(16) 
(17) 



At time t after outburst, the radial and tangential components of 
the magnetic field in a fluid element with density p(R, t) at the 
downstream position R are given by 



B r {R,f) = B d Xti) 



R(t) 



-2 



B t (R,t) = BdM 



p(R,t) R(t) 
Pd(td R s (td ' 



(18) 
(19) 



where f, is the earlier time when this fluid element was shocked 
and pd(ti) is the immediate postshock density at that time. The 
total magnetic field is then simply 



B(R,t) = (B r (R,t) 2 + B,(R,t) 2 ) 



1/2 



(20) 



The radial profile of the total advected magnetic field is 
shown in Figure [6J3. Under the assumption of self-similarity, the 
plotted ratio B(R)/Bj is independent of time. 

2.5.2. Damped downstream magnetic field 



Pohl et al. (2005) suggested that the nonthermal X-ray filaments 
observed in Galactic SNRs could be localized enhancements of 
the magnetic field in the blast wave region. In this scenario, 
the turbulent magnetic field amplified in the shock presursor is 
thought to be rapidly damped behind the shock front by cascad- 
ing of wave energy to very small scales where it is ultimately 
dissipated. 

Assuming a Kolmogorov-type energy cascade of incom- 
pressible magnetohydrodynamic (MHD) turbulence, the charac- 
teristic damping length can be estimated to be (Pohl et al. 2005; 
see also Cassam-Chenai et al . 120071 1 



7 ~ 5 Ud J 

'damp ~ ^max 

n c A ,d 



5 M AM 
0.837T r V2 



(21) 



where uj = V s /r tot is the downstream flow speed in the shock 

1/2 

rest frame, c A ,d = 0.83r tot ca,« is the immediate postshock 
Alfven speed, and A max is the maximum wavelength of the mag- 
netic turbulence, which is expected to be of the order of the gyro- 
radius of the maximum energy protons. In a size-limited shock, 
where the proton maximum energy is determined by upstream 
particle escape, we have 



, _ IfescRjVs 

^max ~ ^V,max — 

?7mfpC 

such that 

'damp 

R s 



3/2 

r tot 7mfp c 



(22) 



(23) 



For V s = 3.35 x 10 4 (f/l day)'"-' km s" 1 (Eq. H, / esc = 0.05, 
T/Wp = 3, and r tot « 4 (Sect. 2.4), one gets 



^~1.3xKT 3 M,„ 



1 day 



We see that for M Am ~ 100, such a magnetic field damping is ex- 
pected to produce thin radio filaments, which were not observed 



in SN 1993J (e.g., Bietenholz et al. 120031 . Therefore, I do not 
use this model of turbulence damping. 

Pohl et al. (2005) also considered that downstream magnetic 
field damping can result from cascading of fast-mode and Alven 
waves in background MHD turbulence. In these cases, we have 
(see also Cassam-Chena'i et al. 2007il 



1 



damp 



2 V27T C A ,d 



V^max^ ■ 



(25) 



where L is the outer scale of the pre-existing MHD turbulence. 
Assuming that this quantity is of the order of the diameter of the 
dense wind bubble blown the red supergiant progenitor of the 
SN, L ~ 2R 0Ut = 6 x 10 17 cm, we obtain 



•damp 

~rT 



0.1M A , U 



1 day 



-1/2 



(26) 



This damping length is thus always larger (for t < 3100 days) 
than the one estimated from the general model of Kolmogorov- 
type energy cascade (Eq.l24l. In the case of damping by cascad- 
ing of fast-mode and Alfven waves the spatial relaxation of the 
postshock magnetic field follows an exponential decay law 



B(R) = B d exp[ - 



R s -R 

^damp 



Postshock magnetic field profiles obtained from Eqs. 
(|27T i with M Am = 9.5 are shown in Figure^. 



(27) 
and 



2.6. Nonthermal electron energy losses 

The cooling processes that can affect the energy distribution of 
the nonthermal electrons during the SN expansion were stud- 
ied by Fransson & Bjornsson dl998l l. They found synchrotron 
cooling to be predominant for electrons radiating at short wave- 
lengths for most of the time, Coulomb cooling to be potentially 
important at early epochs post-outburst, adiabatic cooling to be 
dominant for electrons radiating at 20 cm at late epochs and in- 
verse Compton losses due to electron scattering off photons from 
the SN photosphere to be less important. These results depend, 
however, on several model parameters, e.g. the progenitor mass 
loss rate Mrsg and the magnetic field B u q. 

With the best parameter values of the present model, 
Coulomb cooling is found to be more important than syn- 
chrotron cooling only at very early epochs, when the radio emis- 
sion is still optically thick to SSA (Sect. 4.2). At that time, the 
energy losses of the shock-accelerated electrons have little if no 
effect on the radio emission, such that Coulomb cooling can be 
safely neglected. This contrasts with the model of Fransson & 
Bjornsson (fl998T >. 

I use the work of Reynolds ( 119981 see also Cassam-Chenai 
et al. 120071 Appendix B) to calculate the downstream evolution 
of the electron energy distribution due to synchrotron, inverse 
Compton, and adiabatic losses. The number density of nonther- 
mal electrons per unit energy interval at kinetic energy £ in a 
fluid element being at the downstream position R at time t can 
be written as 



(24) N e {E,R,t) = N e (Ei,R s (ti),td 



p(R,t) 
Pd(td 



4/3 



(28) 



where N e (Ej, R s (ti), f,) is the electron distribution function at the 
initial energy E, produced at the shock at the earlier time f, when 
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Fig. 7. Relative brightness of the 3.6 cm emission from SN 1993 J 
at day 2787 after outburst, as a function of angular radius. In the 
calculations, the outer angular radius 9q is set equal to R s /D. 
Panel (a) shows two unconvolved profiles calculated for / rev = 
(i.e. no contribution from electrons accelerated at the reverse 
shock) and / rev = 30%, assuming in both cases that the magnetic 
field is advected behind the blast wave, with Z?„o = 50 G (B u = 
18 mG at day 2787) and p abs = 0.45R S . Panel (b) shows the two 
profiles convolved with the instrumental resolution of 0.70 mas 
FWHM (0.15 0q at da y 2787) compared to the VLBI data of 
Bietenholz et al. J2Q03I ). 

the fluid element was shockecQ. The change in energy during the 
expansion is given by 



1/3 



Ei 



1 + 0£, 



(29) 



Here, is a radiative loss term that includes both synchrotron 
and inverse Compton cooling: 



0=1^ 



3 (m e c 2 ) 2 



BHt) 
8tt 



+ U md (T) 



qt^ 3 (t)c/t , 



(30) 



where <x x is the Thomson cross section, B(t) is the time- 
dependent magnetic field in the fluid element, and 



U md (T) 



^bol(T) 

4ncR 2 (r) 



is the energy density in the radiation field, Lboi being the bolo- 
metric luminosity. 

Fransson & Bjornsson (1998) showed that during the first 
~100 days after outburst, the radiation energy density at the 



2 We note that in Reynolds ( 1998 1 a factor 07 = p(R, is miss- 

ing in the unnumbered equation just before eq. (24) and in eq. (25). 
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Fig. 8. Calculated radial profiles of the relative emissivity at 
3.6 cm and (f/1 day) = 2787, for B u0 = 30 G and / rev = 10% 
(green dashed lines), B u q — 50 G and / rev — 30% (blue solid 
lines), and B u q = 80 G and / lev = 40% (red dotted lines). The 
shell of uniform emissivity near the contact discontinuity is in- 
tended to account for the emission produced by electrons accel- 
erated at the reverse shock (see text). 



forward shock position was dominated by emission from the 
SN ejecta. They approximated the bolometric luminosity of 
SN 1993 J at these early epochs by 



^bo 



4x 10 



,42 



10 days 



-0.9 



erg s 



(32) 



The ratio of the energy densities of magnetic to seed photon 
fields in the immediate postshock region is then 



Ub 



87Tf/ ra 



B 



i/O 



40 G / \ 1 day 



0.56 



(33) 



This result shows that Compton cooling can play a role soon 
after outburst depending on the magnetic field strength, but that 
synchrotron cooling is expected to become more important after 
some time. A similar result was found by Chevalier & Fransson 
(120061) for Type Ib/c SNe. 

2. 7. Radio synchrotron emission 

Once the electron energy distribution and the magnetic field as 
a function of downstream position are determined as explained 
above, the synchrotron brightness profile and the total flux den- 
sity at a given time can be calculated as described in Appendix B. 



(31) 3. Results 



For a given assumption about the evolution of the postshock 
magnetic turbulence (advection or damping), the model has four 
main parameters to be fitted to the radio data: the mass loss 
rate of the red supergiant progenitor, Mrsg, the normalization 
of the upstream magnetic field, B u q, and the proton and elec- 
tron injection parameters, n p . and rf . . In addition, the radius 

J v ' 'inj 'inj ' 

of the absorbing disk, p a b s , that accounts for FFA of the radio 
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Fig. 9. Same as Fig. 7 but for B M = 30 G, / rev = 10% (green 
dashed lines) and B„o = 80 G, / rev = 40% {red dotted lines). 



emission by the inner SN ejecta (see Appendix B) mainly in- 
fluences the brightness profile as a function of angular radius. 
We first compare calculated synchrotron profiles with the high- 
resolution profile at 3.6 cm measured by Bietenholz et al. (2003 ) 
from VLBI observations. We then study the radio light curves 
reported by Weiler et al. (120071) . 



3. 1 . The radial brightness profile 

To study the radial brightness profile of SN 1993 J with the high- 
est angular resolution, Bietenholz et al. (2003 ) produced a com- 
posite image at 8.4 GHz from VLBI observations performed at 
t = 2080, 2525, and 2787 days after explosion. The data were 
appropriately scaled to take into account the SN expansion and 
then averaged. The resulting brightness profile versus angular ra- 
dius is shown in Figure [7] together with two calculated profiles. 
Here and in the following, we set the radius of the inner opaque 
disk accounting for FFA of the radio emission from the side of 
the shell moving away from us (Appendix B) to p a b s = 0.45/?.,. 
We see in Figure [7J that this simple model of absorption by the 
SN ejecta is consistent with the data for 6/6q > 0.4, but un- 
derestimates the observed emission for 6/Oo < 0.4. The differ- 
ence might arise from incomplete absorption of the radio waves 
in the inner ejecta, possibly because the latter are filamentary 
(see Bietenholz et al. 12003b . However, the excess of emission at 
6/ Go ~ 0.2 accounts for only ~3% of the total flux density and it 
has been neglected so as to limit the number of free parameters. 

The radial brightness profile is mainly sensitive to the 
strength and profile of the magnetic field behind the blast wavfl 
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Fig. 10. Same as Fig. 7 but for the damping of the magnetic tur- 
bulence downstream (green dashed lines) compared to the ad- 
vection of the shocked plasma (blue solid lines). In both cases, 
B M = 50 G and / rev =30%. 

The two calculated profiles shown in Figure|7jare for B u q = 50 G 
and the case of pure advection of the field in the downstream 
plasma. We see in Figure [7p that the model with only electrons 
accelerated at the forward shock (red dashed curve) produces a 
peak at 6 /do ~ 0.85 that is both slightly shifted and too narrow 
as compared to the observed one. There is clearly a deficit of 
emission at 9/Gq < 0.8, i.e. from a region close to the contact 
discontinuity that marks the border between shocked ejecta and 
shocked CSM (see Fig.0. 

As shown in Figure [8] the radial emissivity of electrons ac- 
celerated at the blast wave is cut off before it reaches the con- 
tact discontinuity. This is due to the strong radiative losses suf- 
fered by the electrons that have been accelerated at the earliest 
epochs. The cutoff exists whatever the magnetic lied strength 
B„o; it is mainly produced by Compton cooling for B u o i 10 G 
(see Eq.[33l. Thus, the missing component of synchrotron emis- 
sion is probably coming from another source of accelerated elec- 
trons, most likely the reverse shock. The efficiency of particle ac- 
celeration at reverse shocks in SNRs is poorly known, because 
it depends on the unknown amplification of the magnetic field 
in the unshocked ejecta material (Ellison et al. 2005 ). However, 
recent observations of SNRs provide clear evidence for syn- 
chrotron emission associated with reverse shocks (see Helder & 
Vink 120081 for Cassiopeia A). Here, we simply model this emis- 
sion by a shell of uniform emissivity between the reverse shock 



3 In contrast, the radio morphology of older SNRs, such as Tycho's 

1 



remnant, also depends on the proton injection parameter, T]f_ : , because 



in these objects an increase in the injection efficiency produces a nar- 
rower interaction region (Cassam-Chena'i et al. 20071. This is not the 
case in SN 1993J because of the much higher magnetic field strength in 
the shock precursor region, which limits the backreaction of the shock- 
accelerated nuclei on the shock structure (see Sect. 2.4). 
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and the contact discontinuity (Fig. [8}, whose normalization is fit- 
ted to the radial brightness profile measured by Bietenholz et al. 
(2003 ). We use as normalization factor the integral of the radial 
emissivity relative to that for the forward shock: 



/re 



(Rcd - Rrs)^ S 



Rcd 1 



(R)dR 



(34) 



where Rqd, Rrs, an d ^fs are the radii of the contact discontinu- 
ity, reverse and forward shocks, respectively. We see in Figure|7] 
that for 9/6q > 0.4 an excellent fit to the measured brightness 
profile is obtained for B„q — 50 G and / rev = 30%. After in- 
tegration along the line-of-sight, the synchrotron emission from 
electrons accelerated at the reverse shock is found to contribute 
for 17% of the total detected flux density. 

Figure [8] also shows that the radial emissivity profile de- 
creases more rapidly behind the blast wave as B u q increases, 
which is due to the synchrotron energy losses. Brightness pro- 
files calculated for B u q — 30 and 80 G are compared to the data 
in Figure[9] We adjusted the value of / rev so as to match the maxi- 
mum of the relative brightness to the observed one at Q/Gq = 0.8. 
But we see that the width of the broad emission peak is not 
well fitted: the calculated emission profile is too narrow (resp. 
too broad) for the low (resp. high) value of B u q. Thus, the high- 
resolution brightness profile of Bietenholz et al. provides a first 
indication of the magnetic field strength in the blast wave region: 
1 1 < B u < 29 mG at day 2787 post-outburst (30 < B M < 80 G). 

The measured profile also constraints the evolution of the 
postshock magnetic field. Figure [TUl shows that, when the mag- 
netic field is damped behind the blast wave, the brightness angu- 
lar distribution is shifted toward the outer edge of the radio shell 
and the resulting profile is clearly not consistent with the data. 
This conclusion is independent of both the value of B u and the 
reverse shock contribution. In particular, the brightness profile 
becomes clearly too broad when / rev is increased so as to set the 
maximum of the emission at 6/6q - 0.8. 

The contact interface between the SN ejecta and the shocked 
CSM is thought to be Rayleigh-Taylor unstable and it has been 
suggested that the associated turbulence can amplify the post- 
shock magnetic field (Chevalier et al. 119921 Jun & Norman 
1996). This would increase the synchrotron losses near the con- 
tact discontinuity. Consequently, the position of the cutoff in the 
radial emissivity would likely be shifted to larger radii with re- 
spect to the case with no amplification of the postshock magnetic 
field (Fig. [8}. As a result, a larger contribution of synchrotron ra- 
diation from electrons accelerated at the reverse shock would 
probably be needed to reproduce the measured brightness pro- 
file. 

However, the shocked CSM is expected to be strongly mag- 
netized due to the field amplification by the cosmic -ray stream- 
ing instability operating in the forward shock precursor and the 
postshock magnetic field profile is expected to be dominated by 
the tangential component (see Eqs. 18 and [T9l. Such a strong 
magnetic field could decrease the growth of the Rayleigh-Taylor 
instability and limit a possible additional amplification of the 
field by the turbulence associated with this instability (Jun et 
al. l 19951 >. This aspect of the magnetohydrodynamic evolution of 
the postshock plasma certainly deserves further studies. Here, 
we assume that the postshock magnetic field is not significantly 
amplified by the Rayleigh-Taylor instability. 



3.2. Radio light curves 

In the modeling of the radio light curves, we take into account 
the emission from electrons accelerated at the reverse shock as- 
suming / re v = 30% at all times. As discussed above, this emis- 
sion is then estimated to contribute to a maximum of 17% of 
the total flux density. This number is an upper limit, because the 
synchrotron radiation from the inner shock is strongly attenuated 
when the emission associated with the forward shock is optically 
thick to SSA. We also assume p ar)S = Q.45R S at all epochs. The 
resulting uncertainty on the flux density is negligible. 

Light curves calculated with the best-fit model are shown in 
Figure QT] They were obtained forM RS G = 3.8 x 10~ 5 M Q yr -1 , 
B M = 50 G, 77?. = 10~ 4 , rf^ = 1.1 x 10~ 5 , and pure advection of 
the postshock magnetic field. We see that the general features of 
the data are fairly well represented. But the associated reduced 
X 2 is^fy = 7.9 (x 1 = 4330 for 554 degrees of freedom), which 
is comparable to the one of the best-fit semi-empirical model of 
Weiler et al. ©07): xl = 8.fl However, in the comparison of 
the x 2 values, one could note that the formalism of Weiler et al. 
has nine free parameters (Appendix A), against only four in the 
present model. Interestingly, the present model is also good at 
representing the 12 data points at 0.3 cm (y 2 = 15.5), which is 
not the case for the parametric model of Weiler et al. (x 2 = 159; 
see Fig. [2j). The 0.3 cm flux density between ~10 and ~40 days 
after outburst is lower in the present calculations than in previ- 
ous ones, because the strong synchrotron losses suffered by the 
radiating electrons during the early SN expansion are now taken 
into account. The difference is more pronounced at the short- 
est wavelength, because this emission is produced by higher- 
energy electrons (on average), for which the synchrotron cool- 
ing time is lower (r syn oc 1 /T, F being the electron Lorentz fac- 
tor; see Sect. 4.2). Obviously, this effect cannot be accounted for 
by models neglecting the electron energy losses and assuming a 
homogeneous shell of emission (Chevalier 1 1 982bl 1 1 99 8 ; Perez- 
Torres et al. |200n Weiler et al. 120021 Soderberg et al. [2003T ). 

However, significant deviations of the best-fit model from 
the data can be observed in Figure QT| In particular, we see 
that the straight rising branches of the calculated light curves 
do not represent well the observations at 2 and 3.6 cm. This 
is most likely due to our treatment of external FFA (Sect. 2.2). 
It is possible that the structure of the CSM was more compli- 
cated at the time of explosion than that implicitly assumed by 

clumps cluniDS 

adopting an attenuation of the form (1 - e~ T csM )/r CSM with 
6' = -3m. This would be the case if, for example, the filling fac- 
tor of clumpy material was not constant throughout the whole 
CSM. It is also possible that the CSM temperature 7csm was not 
uniform but varied with radius. However, we note that the data 
at 6 and 20 cm are very well fitted. The best-fit normalization 
to the FFA optical depth is Kj, = 3.0 x 10 4 , which gives from 
Eq. © M RSG = 3.8 x 10~ 5 M yr 1 . 

Like for the brightness profile, the light curves provide clear 
evidence that the postshock magnetic field is essentially ad- 
vected behind the shock. Indeed, we see in Figure QT| that the 
synchrotron flux declines much too rapidly after about day 100 
in the model with damping of magnetic turbulence. 

The modeling of external FFA being uncertain, we now focus 
on the optically thin parts of the radio emission. In Figures Q~2] 
andQj] the flux density is multiplied by a time-dependent power 
law to set upright the decreasing parts of the light curves and the 



4 These high values of x\ are partly due to a few data points that 
strongly deviate from neighbouring points and appear to be incompati- 
ble with the general trend observed at each wavelength. 
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Fig. 11. Radio light curves for SN 1993J at 0.3, 1.2, 2, 3.6, 6, and 20 cm. The data are from Weiler et al. (120071 1 and references 
therein. The blue solid lines represent the best-fit model, which is obtained for Mrsg = 3.8 x 10~ 5 M Q yr , B u q = 50 G, 77^ = 10~ 4 , 

rf m - = 1.1 X 10~ 5 , and pure advection of the postshock magnetic field. The green dashed lines show the case of magnetic turbulence 
damping in the downstream plasma. 



vertical scale is expanded. Figure Q~2] shows the effect of chang- 
ing the magnetic field strength. For each value of B„q, the elec- 
tron injection rate was adjusted to provide a decent fit to the data 
at 3.6 cm. We see that only the model with B u q — 50 G repre- 
sents the data at all the other wavelengths reasonably well. For 
example, for B„o = 100 G the calculated flux densities at 0.3 cm 
fall short of the data, whereas those at 20 cm are too high in 
the optically thin phase. This effect is due to the process of syn- 
chrotron cooling, whose rate increases with the magnetic field 
(T syn oc l/B 2 ) and which steepens the energy distribution of the 
accelerated electrons during their advection downstream. Thus, 
the mean propagated spectrum of nonthermal electrons is too 
steep (resp. too hard) for B„q = 100 G (resp. B„q = 20 G) to 
provide a good fit to the data at all frequencies. It is remarkable 
that when the synchrotron losses are taken into account, the de- 
generacy for the optically-thin emission between the magnetic 
field and the nonthermal electron density (e.g. Chevalier 1998 ) 
is lifted. 

Figure [13] shows the effect of changing 77?.. The values of 
the electron injection parameter are somewhat arbitrary, except 
for the case 77?. = 10~ 4 . We see that in the test-particle case 

= 10~ 5 ), the decline of the optically-thin emis- 



sion with time is too slow as compared to the data. This pro- 
vides evidence that the pressure from accelerated ions is im- 
portant to the structure of the blast wave. But we also see that 
for 7/f j = 3 x 10~ 4 the light curves decline too rapidly, which 
shows that the shock modification is relatively weak. There is 



a partial correlation between 77^ and B u q. By varying these two 
parameters and comparing the calculated curves to the data, I 
estimate the acceptable range for the proton injection rate to be 
5 x 10~ 5 < 77?. < 2 x 10~ 4 . Higher values of 77?, are often used 
in studies of Galactic SNRs (e.g., Cassam-Chenai et al. (2007) 
adopted 77^ = 10~ 3 for Tycho's remnant, based on the observed 
closeness of the blast wave and contact discontinuity). 

The time evolution of the synchrotron flux is related to 
the energy distribution of the emitting electrons (see, e.g., 
Chevalier 1982b): the steeper the electron spectrum, the faster 
the decline of the optic ally- thin emission. As dicussed in 
Sect. 2.4, the main effect of the cosmic-ray pressure is to reduce 
the compression ratio of the subshock, r su b, whereas the overall 
compression ratio r m remains nearly constant. This shock mod- 
ification affects essentially the particles of energies < m p c 2 that 
remain in the vicinity of the subshock during the DSA process. 
As shown in Figure [4] an increase of 77^ causes the energy dis- 
tribution of the nonthermal electrons to steepens below ~ 1 GeV, 
which is the relevant energy domain for the radio synchrotron 
emission. The light curves provide evidence that the spectral 
index of the shocked electrons in this energy range is slightly 
higher than the test-particle value q su \, = 4, which is effectively 
obtained with 77 P . * 10 -4 . 

'inj 

Although in most previous models of radio emission from 
SNe the shocked-electron spectral index is assumed to remain 
constant, a gradual steepening with time of the shocked-particle 
energy spectrum below m p c 2 is to be expected. This is because of 



14 



V. TatischefF: Radio emission and nonlinear diffusive shock acceleration in SN 1993J 




Days after outburst Days after outburst 

Fig. 12. Radio light curves of SN 1993 J with the flux density multiplied by (f/f^) 7 . The values of tx for each wavelength are 
indicated in the Figure. Green dashed lines: B u q — 20 G and 77^. = 2.3 x 10~ 5 ; blue solid lines: B M — 50 G and rf^. — 1.1 x 10~ 5 ; 

red dotted lines: B u0 = 100 G and 77 s . = 8 x 10~ 6 . The other parameters are M RSG = 3.8 x 10~ 5 M yr _1 and if . = 10~ 4 . 



the decrease of the upstream magnetic field, which reduces the 
efficiency of gas heating via Alfven wave damping in the precur- 
sor region and consequently reduces the subshock compression 
ratio. Thus, for rf m - = 10~ 4 , we find that r su b decreases from 3.89 
to 3.70 between day 10 and day 3100 after outburst (Fig. [5]), 
which causes q m b to increase from 4.04 to 4.1 1 (Eq.ll4l. 

For given values of B u q and 77?., the normalization of the 
optically-thin synchrotron flux is set by the product p u x if - °c 
Mrsg x rf m y As the progenitor mass loss rate also fixes the 
level of FFA in the CSM, the electron injection rate can be 
uniquely determined. For 77?. = 10~ 4 , the best-fit parameter is 

77^. = 1.1 x 10~ 5 . To set the amplitude of the nonthermal elec- 
tron spectrum in DSA models, several authors use the electron- 
to-proton density ratio at relativistic energies, K ep , rather than the 
electron injection rate (e.g., Ellison et al. 2000). Here, we obtain 
at, e.g., 10 GeV K ep = 2.2 x 10~ 3 at day 10 and K ep = 1.8 x 10~ 3 
at day 3100. The value of K ep is higher at day 10 because g sl ,b is 
lower (see Eq.[T3l. The obtained values of K ep are roughly con- 
sistent with those estimated for the blast wave of Tycho's SNR: 
4xKT 3 (Volk et al. 120021 -1-5 x 10~ 3 (Cassam-Chenai et al. 
[20071 . 

4. Discussion 

The model gives a consistent description of the large set of 
radio data available for SN 1993 J. In particular, it is remark- 
able that the same magnetic field strength, B„q ~ 50 G (in 
B u = B u o{t/l day) -1 ), is required to explain the high-resolution 



brightness profile measured by Bietenholz et al. (120031 ) at day 
2787 after outburst and the set of radio light curves reported by 
Weiler et al. (120071 ). This further supports the assumption that the 
magnetic field at the shock varies approximately as l/t. Besides, 
both the brightness profile and the light curves show that the 
postshock magnetic field is not damped behind the blast wave 
but simply carried along by the plasma flow in the downstream 
region. We shall see in Sect. 4.3 that the obtained magnetic field 
evolution can be explained in terms of an amplification of the 
fluctuating component of the field by the Bell (2004)'s nonreso- 
nant cosmic-ray streaming instability in the shock precursor re- 
gion. But before, we are going to discuss in the light of the re- 
sults two important assumptions of the model: the adopted CSM 
density profile (Sect. 4.1) and the relevant processes of cooling 
of the accelerated electrons (Sect. 4.2). 

4.1. Circumstellar density profile and mass loss rate of the 
progenitor star 

As already discussed in Sect. 2.2, there has been considerable 
debate about the stru cture o f the CSM of SN 1993 J. Several pa- 
pers (Van Dyk et al. 1994; Fransson et al. 1996; Perez-Torres 
et al. I200T1 Immler et al. I200T1 Weiler et al. 12007] ) advocated 
for a flatter density profile, pcsm 00 R " with s = 1.5-1.7, than 
the s — 2 power-law expected for a constant mass-loss rate and 
constant-velocity stellar wind prior to the SN explosion. On the 
other hand, Fransson & Bjornsson d 19981 ) obtained a good de- 
scription of the radio data using the standard s - 2 density pro- 
file. 
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Fig. 13. Same as Fig. 12 but for various values of the injection parameters: 77?. = tjf^ = 10 5 {red dotted lines); . = 10 4 and 
= 1.1 x 10~ 5 {blue solid lines), and 17?. = 3 x 10~ 4 and 77^ = 1.8 x 10~ 5 {green dashed lines). The other parameters are 
Mrsg = 3.8 x 10~ 5 M yr 1 and B M = 50 G. 



Figure [14] shows light curves calculated with s = 1.6. The 
two sets of theoretical curves represent two different assump- 
tions for the postshock magnetic field evolution: Bd cc £ _1 , con- 
sistent with the analysis of Sect. 2.3, and Bd cc f~ . The latter 
case, which is shown for completeness, can result from either 
a scaling of the magnetic energy density with the total post- 
shock energy density {B 2 d cc p u V 2 ; like in models 1 and 3 of 
Chevalier 1996) or a scaling of the postshock magnetic field 
with the pre-existing circumstellar field {Bd cc R l ; like in mod- 
els 2 and 4 of Chevalier 1 19961 1 . In order to reproduce the ob- 
served rate of decline of the radio emission in the optically thin 
phase, the spectrum of the emitting electrons has to be steeper 
for the flatter CSM density profile than for s - 2. This is be- 
cause F v cc N e cc p u . It implies that rf. should be higher for 
s = 1.6 than for s = 2. But we see that the resulting distri- 
butions of shocked electrons cannot provide simultaneously a 
good fit to the data at all wavelengths. The theoretical flux den- 
sities are clearly too low at 0.3 and 1.2 cm and too high at 6 
and 20 cm as compared to the observations in the optically thin 
phase. This is somewhat similar to the case of the too high mag- 
netic field shown in Figure Q~2] But with the low field strength 
B u o = 10 G adopted here, synchrotron losses are less important 
than Compton losses at early times (see Eq. l33l and than adia- 
batic losses latter on. Therefore, reducing B u q below 10 G (and 
increasing rjf^ at the same time ) does not help to improve the fit. 
Thus, we conclude that it is not possible to reproduce the obser- 
vations with s = 1.6 in the framework of the present model. 

Weiler et al. (2007 ) recently adopted s - 1.61 based on pre- 
vious works and interpreted the flatter CSM density profile in 



terms of a steady decrease of the mass loss rate of the progenitor 
star, from 5.9 x 10~ 6 M Q yr ~8000 years before explosion to 
5.4 X 10" 7 M G yr at the time of the SN. These values are much 
lower than the mass loss rates estimated in the present as well as 
in previous works, which are in the range (2-6)xl0~ 5 M yr~' 
just prior to explosion (Van Dyk et al. |1994| Fransson et al.] 1996 ; 
Fransson & Bj6rnsson |19981 Immler et al. 120011 ). 

As discussed in Appendix A, the use of the fitting formalism 
developed by Weiler et al. (Q986 2002 and references therein) 
may lead to uncertain physical parameters when both SSA and 
FFA in the CSM are important to explain the radio light curves. 
In this case, it could be justified to fix in the fitting proce- 
dure the time dependence of the mean magnetic field in the 
synchrotron-emitting region. As a check, I performed a simul- 
taneous, least-squares fit to the light curve data for SN 1993 J 
using the Weiler et al.'s formalism together with Eq. (IA.13I ). that 
is fixing 5" — B - 2m - 0.5. 1 found the attenuation of the radio 
flux by the homogeneous component of the CSM to be not sig- 
nificant (i.e. K2 compatible with zero) so I also fixed A^™ 8 = 1 . 
I then obtained the best-fit parameters K\ = 5.4xl0 3 , a = -0.79, 
B = -0.75, K 3 = 1.16 x 10 5 , 5' = -2.42, and K 5 = 6.9 X 10 4 , 
with x 2 = 4250 (for 552 degrees of freedom). The value of K3 
gives from Eq. ([8]( Mrsg = 7.6 x 10~ 5 M Q yr~', which is within a 
factor of two of the result obtained from the detailed model pre- 
sented above. Noteworthy, the best-fit value of 6' is very close to 
the theoretical expectation: 6' — m{\ — 2s) = -2.49 for s = 2. 
This provides support to the assumption that the external free- 
free attenuation of the radio emission is mainly due to the inho- 
mogeneous component of the CSM. 
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Fig. 14. Same as Fig. 1 1 but for a CSM density profile of the form pcsm <* R 1 6 - Green dashed lines: B u oc t 1 
red dotted lines: B u oc f-°- 83 and rf m . — 5 x 10 -3 . For simplicity, FFA in the CSM is taken into account as before, with 5' = -3m and 
K 3 = 3.0 x 10 4 . The other parameters are B uQ = 10 G and rf. = 10" 2 . 



4.2. On the Coulomb energy losses of the nonthermal 
electrons 

The coulomb energy loss rate of relativistic electrons is given by 
(see Gould [19751 
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where F cou \ = ln(T ' 2 /e p ) + 0.216, cr x is the Thompson cross 
section, n e d is the number density of thermal electrons in the im- 
mediate postshock plasma, 
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and e p = 7.27 X 10~ u n e d (in cgs units) is the normalized plasma 
energy. Coulomb cooling is more efficient than synchrotron 
cooling when r cou i < T syn , the synchrotron energy loss rate being 
given by 
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The typical Lorentz factor for the relativistic electrons radi- 
ating at frequency v can be written as (see Eq. 
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where we have used for the magnetic field configuration across 
the shock B4 * 0.83r tot B u (see Eqs. 16 andfTTli. Inserting F c into 
Eq. ( f35t and ( [37] >, one obtains a limit on the time during which 
Coulomb cooling is more important than synchrotron cooling: 
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Here, we have used F cou i ss 30 and X = 0.3. This result shows 
that Coulomb energy losses can be safely neglected for frequen- 
cies v > 8.4 GHz. They could play a role, however, during the 
rising phase of the light curve at 1 .4 GHz (20 cm), for ~ 210 days 
after outburst. But at that time, the 20 cm emission was attenu- 
ated by SSA and the emitted flux density depended only weakly 
on the propagated spectrum of nonthermal electrons. Thus, the 
neglect of Coulomb cooling in the model is justified. 

4.3. Magnetic field amplification 

Fransson & Bjornsson (1998) modeled the radio emission of 
SN 1993 J assuming that the radiation is emitted from a shell of 
uniform magnetic field. They found {B) w 340(f/l day) -1 G (or 
alternatively (B) * 64(R S /10 15 cm) -1 G), which is in reasonable 
agreement with both our first estimate (Eq. |T0T > and the mean 
postshock magnetic field resulting from our detailed analysis. 
Our best-fit value for the immediate postshock magnetic field is 
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B d * 0.83r tot B„ * 180(f/l day)-' G, given r tot * 4.3 (Fig.0. 
The mean magnetic field strength in the synchrotron-emitting re- 
gion is expected to be higher than Bj by a factor of less than two, 
given that the field is mainly advected in the downstream plasma 
(see Fig. 01, i.e. 180 < (B)(t/1 day) < 360 G. We note, however, 
that in our detailed model all the emitting electrons do not see 
the same mean magnetic field, because due to the energy losses 
the downstream transport of the nonthermal particles depends on 
the energy they acquired at the shock front. 

Perez-Torres et al. (1200 11 1 were able to obtain a good fit to 
the radio light curves with {B) * 56(?/l day) -0 86 G. But they 
did not take into account the electron energy losses and assumed 
the existence of a low-energy cutoff in the distribution of the rel- 
ativistic particles, whose main effect is to flatten the radio light 
curves at the shortest wavelengths. But as discussed in Sect. 3.2, 
the observed flattening of the light curves at early epochs is nat- 
urally explained by the strong synchrotron losses suffered by the 
radiating electrons. Moreover, a low-energy cutoff in the particle 
distribution function is not supported by the DSA theory. 

As already pointed out by Fransson & Bjornsson ( 1998) and 
Perez-Torres et al. (2001) the magnetic field inferred from the ra- 
dio observations is at least two orders of magnitude higher than 
the magnetic field of stellar origin that is expected to preexist 
in the wind of the red supergiant progenitor. The high magnetic 
field in the blast wave region is likely a result of an amplifi- 
cation associated with the production of nonthermal particles 
by nonlinear DSA. The diffusive streaming of cosmic-rays in 
the upstream plasma is expected to produce strong turbulence, 
which could amplify the chaotic component of the magnetic field 
(e.g., Bell & Lucek I200T1 Amato & Blasi 120061 Vladimorov 
et al. |20061 l . Two different mechanisms of turbulence genera- 
tion are discussed for SNR shocks. First, due to the cosmic -ray 
anisotropy in the shock precursor, a resonant instability is pre- 
dicted to exist in this region and to produce a rapid and intense 
excitation of Alfven waves (e.g., McKenzie & Volk 1982; Bell & 
Lucek 2001} • Bell (|2004) also found a nonresonant streaming in- 
stability caused by the Lorentz force associated with the cosmic- 
ray electric current, which could strongly amplified MHD per- 
turbations of short wavelengths. The overall magnetic field am- 
plification could be the result of both instabilities operating in 
the shock precursor. 

The saturation value of the random magnetic field amplified 
by the nonresonant instability, SB m , can be estimated from equa- 
tion (28) of Pelletier et al. d2U06l see also Bell |2004l l: 
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Fig. 15. (a) Normalized immediate postshock pressures in ther- 
mal (P t h) and nonthermal (Pgr) particles (f tot = P± + Pgr) a s 
a function of time after outburst, for the best obtained solution 
(M RSG = 3.8 x 10- 5 M e yr\ B M = 50 G, ijg. = 10" 4 , and 

1 . 1 x 10~ 5 ). (b) Magnetic field immediately upstream from 
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the subshock B u times (f/1 day). The hatched area shows the 
best-fit value obtained from the radio data, B u = Z?„o(f/l day) -1 
with B u o = 50 + 20 G, and the red solid curve shows the fluc- 
tuating magnetic field expected from the nonresonant instability 
(eq.EOll. 
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where <t> = \n{p max lm p c) ~ 15 (see Sect. 4.4 below) and Pgr 
is the cosmic-ray pressure at the shock front. The normalized 
pressure £cr = Pcr/PuV 2 calculated with the nonlinear DSA 
model from the best-fit parameters is shown in Figure [T5K We 
see that £cr increases with time, but is always much lower than 
the normalized postshock thermal pressure. For such a weakly 
cosmic-ray-modified shock, we have in good approximation 
£ C r <* p'yv 2 s (see Berezhko & Ellison [TP991 eq. 34), which 

leads to £cr k V' 1 given that . oc p th oc V s (Blasi et al. 2005). 
Here, /? t h is the most probable momentum of the thermal protons 
in the postshock plasma. We then obtain from Eq.d40li 



by Bell (2004) is consistent with that inferred from the data. As 
shown in Figure [T5b, the strength of the field amplification cal- 
culated from Eq. (l4Qb is also quite close to that required to ex- 
plain the radio flux. 

Following Pelletier et al. (2006) the relative importance of 
the resonant and nonresonant regimes of the streaming instabil- 
ity can be estimated from 
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As previously stated (Sect. 2.3) the magnetic field evolution ex- 
pected from the nonresonant instability mechanism uncovered 



This expression shows that in SN 1993 J the nonresonant mech- 
anism is expected to dominate at early times, but the resonant 
instability becomes more and more important as the shock is 
slowing down. 
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According to Bell (2004 ) the minimum timescale for growth 
of MHD waves driven by the nonresonant instability at a dis- 
tance z upstream from the subshock is given by 



T ar (z) 



®cE p (z) 



(43) 
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where —e is the electronic charge, e nt is the cosmic -ray acceler- 
ation efficiency (see Figure [5]) and E p (z) « 3zeB u V s c^ l j]^ is 
the energy of the accelerated protons whose upstream diffusion 
length is equal to z. Numerically, the MHD growth timescale for 
SN 1993J is 



T m (z) * 3.3 x 10 l \ — 



4> \( % 
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day. 



(44) 

Thus, we have r nr «: t for E p (z) < 10 15 eV. This shows that 
the turbulent magnetic field amplified by the nonresonant insta- 
bility should reach its saturation value over most of the shock 
precursor region soon after the onset of particle acceleration. 

Because the damping of MHD turbulence can efficiently heat 
the thermal gas in the shock precursor region, a high amplifi- 
cation of the fluctuating magnetic field could at the same time 
provide a strong limitation for the modification of the shock 
structure induced by the cosmic -ray pressure. According to, e.g., 
Berezhko & Ellison i 19991 1 the degree of heating of the shock 
precursor via Alfven wave damping is set by the ratio M\ u /Ma, u 
(Ms,u and Ma, u are the upstream sonic and Alfven Mach num- 
bers, respectively; see Sect. 2.4). Identifying B„ and 6B m one 
obtains from Eq. (f4Tb the scaling relation 



Am 



(45) 



which shows that for the same cosmic -ray pressure the shock 
modification should be less pronounced in a high-speed shock 
than in a slower one. 

Bjornsson & Fransson (2004) modeled the radio and X-ray 
emission from SN 2002ap with an accelerated electron energy 
index y « 2, as expected for a test-particle shock. The above 
scaling relation suggests that the absence of clear shock modifi- 
cation in this object is a consequence of the very high velocity 
of the blast wave (V s ~ 7 x 10 4 km s _1 ; Bjornsson & Fransson 
120041 . 

Various studies of the synchrotron emission produced 
in young SNRs led Berezhko & Volk (2006, see also 
Berezhko 2008 ) to propose the following empirical relation be- 
tween the amplified magnetic field energy and the cosmic-ray 
pressure: 
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For SN 1993 J, a better relation would be (see Eq.l40l 
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It is remarkable that this expression could provide a fair esti- 
mate of the magnetic field amplification in the high-speed blast 
wave of SN 1993J (V s /c ~ 0.1), as well as in slower shocks of 
Galactic SNRs, which are at the end of the free expansion phase 
of the post-SN evolution or in the Sedov phase (VJc ~ 0.01). 
In these objects, however, the resonant regime of the cosmic -ray 
streaming instability is expected to dominate the field amplifica- 
tion. 
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Fig. 16. Maximum energy of shock-accelerated protons as a 
function of time after outburst. E^ x and E s ^ x are the maximum 
energies caused by the finite shock age and size, respectively. 
The true maximum proton energy (solid line) is the minimum of 
these two quantities. 



4.4. SN 1993J and the origin of cosmic rays 



As already pointed out by Bell & Lucek (2001) cosmic -rays 
were probably accelerated to very high energies in the blast 
wave of SN 1993 J short after shock breakout. Calculated maxi- 
mum proton energies are shown in Figure [16] The quantity EjJSx 
is the maximum energy caused by the finite age of the shock 
(Baring et al. 1999; Tatischeff & Hernanz [20071 ) and is obtained 
by time integration of the DSA rate from an initial acceleration 
time assumed to be to = 1 day after outburst. The maximum 
energy E^ x is caused by the finite spatial extend of the shock 
and is calculated by equalling the upstream diffusion length of 
the accelerated protons to some fraction / esc of the shock radius. 
Numerically, we find (see, e.g., Tatischeff & Hernanz 2007 ) 
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(48) 



We see in Figure [16] that for / esc = 0.05 the size limitation of 
the shock becomes rapidly more restrictive than the age limita- 
tion. This conclusion is independent of t^. However, the value of 
/e S c, which depends on the upstream turbulence generated by the 
accelerated particle streaming (Marcowith et al. 2006), is uncer- 
tairQ But in any case, the resulting maximum proton energy (i.e. 
the minimum of E^ x and is lower than the previous esti- 
mate of Bell & Lucek ( T200TT ). ~ 3 x 10 17 eV. Moreover, we 
did not take into account the nonlinear effects recently pointed 
out by Ellison & Vladimirov (2008), which could further reduce 



5 A simple estimate of / csc may be derived by assuming that the up- 
stream escape boundary at z = fe X R s corresponds to the distance for 
which the MHD growth timescale T nr (z) = t. We then obtain from 
Eq. g3) 
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For e n , = 0.1, * = 15, r] mfp = 3, m = 0.83, and M A ,„ = 9.5 (see 
above), we get / csc « 0.05. However, the full interpretation of this result 
is beyond the scope of this paper. 
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E'Lax- However, it is likely that SN 1993J has accelerated protons 
above 3 x 10 15 eV, i.e. the energy of the spectral "knee" above 
which the all-particle cosmic -ray spectrum measured near Earth 
shows a significant steepening. This could provide support to 
the scenario first proposed by Volk & Biermann ( 1988) that the 
very-high energy Galactic cosmic-rays are produced by massive 
stars exploding into their former stellar wind. 

In this context, it is instructive to estimate the total energy ac- 
quired by the cosmic -ray particles in SN 1993 J during the early 
stage of interaction between the SN ejecta and the red supergiant 
wind: 



J la 
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7.4 x 10 49 erg 



(50) 



where tj = 3 100 days and e nt is the acceleration efficiency shown 
in Figure [5] It is remarkable that the value obtained for £cr is 
very close to the mean energy per SN required to account for 
the Galactic cosmic -ray luminosity, a; 7.5 x 10 49 erg (Tatischeff 
2008). Before day ~3100, i.e. during the first ~8.5 years after 
explosion, the forward shock processed in the expansion a total 
energy of 4. Ox 10 50 erg and the mean acceleration efficiency was 
<e nt > = 19%. 

Most of the cosmic rays accelerated during this early stage 
were advected downstream from the shock and will remain in- 
side the SNR until it merges with the interstellar medium. These 
particles suffer adiabatic losses and their energy decreases like 
R~ l (e.g. Drury [19961 Thus, they will make a negligible contri- 
bution to the interstellar cosmic -ray population at the end of the 
SNR. But as shown in Figure[l6] the highest-energy cosmic -rays 
are expected to escape continuously the shock system from up- 
stream, due to the finite size of the acceleration region. The value 
obtained for Eqr (Eq. l50i l suggests that these very-high energy 
particles could be important to the origin of interstellar cosmic 
rays above ~ 10 15 eV. 

At day ~3 100 after outburst, the blast wave started to expand 
into a more diluted CSM resulting from a lower mass loss rate of 
the SN progenitor more than ~9000 years before explosion (see 
Weiler et al. 120071 The CSM density structure of SN 1993J can 
be explained by stellar evolution calculations that show that red 
supergiant stars can undergo a stage of intensified mass loss for 
few 10 4 years before SN explosion (Heger et al. 1997). After day 
~3100 post-outburst, an additional number of very-high energy 
cosmic -rays could have escaped the shock region into the inter- 
stellar medium, because their diffusion length t u has increased 
{t u oc B~ l oc p~^ 2 ; see eq. |4~TT>. Significant variations in the mass 
loss activity of red supergiant stars prior to explosion are com- 
monly inferred from radio observations of type II SNe (Weiler 
etal.|2p02J. 

These results suggest that Galactic cosmic rays of energies 
above ~ 10 15 eV could be produced immediately after the explo- 
sion of massive stars (Volk & Biermann |1988l BiermannQ993), 
provided that the acceleration efficiency found for SN 1993 J is 
typical of most SNe and that the mass-loss rate obtained for the 
progenitor of this type lib SN, M RS g = 3.8 x 10~ 5 M Q yr -1 , is 
close to the mean mass-loss rate of red supergiant stars prior to 
explosion. 



4.5. SN 1993J as a gamma-ray source 



Kirk et al. ( I1995I ) pointed out that radio SNe could emit a large 
flux of very-high energy gamma rays arising from the decay of 
neutral pions created in hadronic collisions of cosmic rays with 
target nuclei. A detection of the predicted TeV emission with 



ground-based atmospheric Cherenkov telescopes would provide 
a direct confirmation of efficient acceleration of high-energy pro- 
tons in SN explosions. 

According to Kirk et al. ( 1995), SN 1993 J was a promising 
target for detection, with a predicted flux of gamma-rays above 
1 TeV of F y (> 1 TeV) = 2 x 10~ 12 cnr 2 s" 1 for years after out- 
burst. But their calculations assumed a much flatter CSM den- 
sity profile, pcsm K R~ 15 , than the standard radial dependence 
Pcsm K R 2 found here. 

The gamma-ray flux can be estimated from 
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where q y (> 1 TeV) a 10~ 17 photons s erg~' cm 3 (H-atom) 1 
is the gamma-ray emissivity, which depends on the cosmic-ray 
spectrum produced at the shock (Drury et al. |1994| l and (n p ) ~ 
n e d /(l+2X) (see eq.l36ll is the mean density of hydrogen nuclei in 
the downstream volume V occupied by the cosmic rays, whose 
mean energy density is (ccr) ~ 3 Pgr- in first approximation we 
can identify the volume V with that of the shocked CSM: 



(52) 



where Rcd = /f. v /l-27 is the radius of the contact discontinuity 
(Fig. [6}. We then obtain 
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(53) 



where we have used £cr ~ 0.04(f/l day) 1- ™ (see Fig. [TBI 

But the TeV gamma-ray emission can be attenuated by the 
pair production process y + y — > e + + e~ in the strong radiation 
field from the SN ejecta. The corresponding optical depth at the 
gamma-ray energy E y can be estimated from 



where 
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is the absorption coefficient in the diluted blackbody radiation 
field of temperature T b b and energy density at the shock £/ ra d 
(Gould & Schreder 1967). The SN ejecta emitted mainly in the 
UV band during the first week after explosion and then in the 
optical for more than 120 days with a mean blackbody temper- 
ature r bb a 7000 K (Richmond et al. [19941 Lewis et al. [19941 
For this temperature, the function f yy (E y , 7bb) reaches its max- 
imum value » 1 at E y m 1 TeV (see Gould & Schreder |19671 l. 
Using t/, a( j « Lboi/(4^c/? 2 ) with the power-law approximation of 
Fransson & Bjornsson (1998) for the evolution of the SN bolo- 
metric luminosity (Eq.[32l we obtain 
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1.73 



(56) 



Thus, the early TeV emission of SN 1993 J was strongly attenu- 
ated by photon-photon absorption. Multiplying the unabsorbed 
flux F 7 (> 1 TeV) given by Eq. ((53) by the attenuation factor 
A yy - exp(-T rr ), we find that the flux reached a maximum of 
»4 x 10~ 15 cm~ 2 s~' about 270 days after explosion, which is a 



20 



V. TatischefF: Radio emission and nonlinear diffusive shock acceleration in SN 1993J 



factor of »40 lower than the 5<x sensitivity of the High Energy 
Stereoscopic System (H.E.S.S.) of atmospheric Cherenkov tele- 
scopes for 25 hours of observation near zenith (Aharonian et al. 
I20CJ6I 

Below ~50 GeV, gamma rays are not attenuated by pair pro- 
duction in collisions with optical photons. The total gamma-ray 
flux in the energy range 1-50 GeV can be readily estimated 
from Eq. ( T5TT > using the emissivity q 7 (> 1 GeV) * 1CT 14 pho- 
tons s _1 erg -1 cm 3 (H-atom) -1 (Drury et al. 11994) . It gives for 
SN 1993J 
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(57) 



The mean flux of gamma rays above 1 GeV from day 1 to day 
8 after outburst is then »;6xl0~ 10 crrT 2 s , which is *14 times 
lower than the 5<x sensitivity of the LAT instrument on the Fermi 
Gamma-ray Space Telescope for one week of observation in the 
survey modqj This suggests that a SN similar to SN 1993 J in 
terms of progenitor wind properties and cosmic-ray acceleration 
efficiency could be detected in 7r°-decay gamma-rays with Fermi 
out to a maximum distance of ~1 Mpc. 

5. Conclusions 

Inspired by previous works on the morphology of synchrotron 
emission in Galactic SNRs (Cassam-Chenai et al. 12005) . I have 
employed a model to explain the radio emission of SN 1993 J, 
which includes nonlinear effects of the diffusive shock accelera- 
tion mechanism. The model states that cosmic-ray ions are accel- 
erated at the SN blast wave and their pressure modify the shock 
structure. Therefore, the energy spectrum of the accelerated elec- 
trons depends on the rate of injection of shocked protons into the 
acceleration process. The nonthermal electrons suffer both adi- 
abatic and radiative losses during their advection downstream 
the shock, which modify their energy distribution. The hydrody- 
namic evolution of the postshock plasma is calculated by assum- 
ing that the SN expansion is self-similar. The relativistic elec- 
trons emit synchrotron radiation in the ambient magnetic field 
and the corresponding radio emission is worked out from radia- 
tive transfer calculations that include the process of synchrotron 
self-absorption. The magnetic field amplification associated with 
the efficient production of cosmic -rays is not calculated in the 
framework of the model, but instead the strength of the amplified 
field is determined from a fit to the radio data. The model con- 
tains three other parameters to be fitted to the data: the shocked 
proton and electron injection rates into the DSA process and the 
mass loss rate of the SN progenitor. 

The extensive radio observations of SN 1993 J make this ob- 
ject a unique "laboratory" to study cosmic-ray acceleration in 
a SN shock. In particular, the VLBI imaging observations that 
resulted in detailed measurements of the SN expansion are of 
prime interest for the DSA theory. By applying the model to 
both the high-resolution brigthness profile at 3.6 cm measured 
by Biethenholz et al. (2003) and radio light curves at six fre- 
quencies reported by Weiler et al. (12007) , I have obtained the 
following main conclusions: 

1. The CSM density profile is consistent with a constant, 
steady wind of the red supergiant progenitor of the SN for 
~9000 years before explosion. Contrary to previous claims, 
SN 1993 J is not markedly different from the other studied 



radio SNe with regards to the mass loss evolution of the 
progenitor star. The best-fit mass-loss rate is Mrsg = 3.8 X 
10~ 5 M Q yr~' for the assumed stellar wind terminal velocity 
u w = 10 km s _1 and CSM temperature 7csm = 2 x 10 s K. 

2. The observed morphology of the radio emission provides ev- 
idence that part of the synchrotron radiation («17% of the 
total flux density at day 2787 post-outburst) is produced by 
electrons accelerated at the reverse shock. 

3. The best-fit rate of injection of cosmic-ray protons into the 
DSA process at the forward shock is 77?. = 10~ 4 , with an 
estimated error of a factor of two. The calculated fraction of 
the total incoming energy flux converted by the blast wave to 
cosmic-ray energy during the first ~8.5 years after explosion 
is (e n t) = 19%. The shock modification induced by the back- 
pressure from the accelerated ions is increasing with time af- 
ter outburst. However, the cosmic-ray pressure accounts for 
less than 20% of the total postshock pressure at all times, 
such that the shock modification remains relatively weak. 

4. The best-fit cosmic-ray electron injection rate is if. = 1.1 X 

10 and the calculated electron-to-proton density ratio at 



10 GeV isK t 



ep 



2x10 3 . The obtained value of K ep is lower 



than that measured in the Galactic cosmic-rays near Earth, 
ssl0~ 2 . The latter number, however, results from cosmic -ray 
transport in the Galaxy and therefore is probably not indica- 
tive of the relativistic electron-to-proton ratio at the source 
of cosmic -ray acceleration. 
5. The synchrotron energy losses suffered by the radiating elec- 
trons in the postshock magnetic field are important for the 
modeling of both the radio light curves and the morphol- 
ogy of the radio emission. The inferred magnetic field is 
broadly consistent with that expected from an amplification 
in the shock precursor region by the nonresonant regime 
of the cosmic-ray streaming instability. The magnetic field 
immediately upstream from the subshock is found to be 
B u = B u0 (t/l day)'' with b * -1 and B M = 50 ± 20 G. The 
measured field strength is higher by a factor of 2 to 5 than 
the saturated value of the turbulent magnetic field predicted 
from the model of Bell (2004) for the turbulence generation. 
The energy density in the amplified magnetic field is found 
to obey the relation 
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6 As calculated with the LAT source detectability tool at 
http://fermi.gsfc.nasa.gov/ssc/proposals/detectability.html 



6. The turbulent magnetic field amplified in the precursor re- 
gion is not damped behind the shock (as proposed by Pohl et 
al. 2005 for Galactic SNRs) but essentially carried along by 
the plasma flow in the downstream region. 

7. The magnetic field amplification increases the DSA rate, 
thus allowing the rapid acceleration of cosmic -ray protons to 
energies well above 10 15 eV. The proton maximum energy is 
find to be limited by the finite confinement size of the shock, 
which implies that the highest-energy cosmic rays continu- 
ously escaped the acceleration region from upstream during 
the first ~8.5 years after explosion. After that time, when the 
blast wave passed the limit of the dense CSM established by 
a stage of high mass loss from the red supergiant progenitor 
prior to explosion, an additional number of very-high energy 
protons probably escaped into the interstellar medium. The 
results obtained for this SN provide support to the model 
of Volk & Biermann ( [1988) and Biermann ( fl993l that the 
Galactic cosmic-rays above ~10 15 eV are accelerated in the 
explosion of massive stars, during the early stage of interac- 
tion of the SN ejecta with the progenitor wind. 
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8. The early emission from SN 1993 J of TeV gamma-rays pro- 
duced in hadronic collisions of accelerated cosmic -rays with 
ambient material was strongly attenuated by pair production 
in the dense radiation field from the SN ejecta. The flux at 
Earth of photons of energy above 1 TeV reached a maxi- 
mum of only =s4 x 1CT 15 cm -2 s _! at day ~270 after explo- 
sion. Above 1 GeV, the flux is found to be F y (> 1 GeV) « 
2 x 10~ 9 (f/l day) -1 crrT 2 s _1 . This result suggests that type 
II SNe could be detected in 7r°-decay gamma-rays with the 
Fermi Gamma-ray Space Telescope out to a maximum dis- 
tance of only ~1 Mpc. 
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Appendix A: Magnetic field evolution from the 
parametric model of Weiler et al. 

Weiler et al. ( 1986 2002 and references therein) have developed 
a semi-phenomenological model to described the light curves of 
radio SNe. In this model, the flux density at a given frequency v 
and time t after outburst can be expressed as 

- K {ji^{x^ A ™ A ™ A ^ ■ (A1 > 

where 

C = -p(-r h c r M ° 8 ) , (A.2) 
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Here, K\ represents the unabsorbed flux density in mJy, and K2, 
Kj,, and K$ are the optical depths for attenuation by a homo- 
geneous absorbing CSM, a clumpy or filamentary CSM, and 
internal SSA, respectively, at day one after outburst and for 
v = 5 GHz. The parameters 5, 6' ', and 6" describe the time 

dependance of the optical depths 8 , T^g™ ps , and tssa- The 
parameters a and B are the spectral index and rate of decline, 
respectively, in the optically thin phase. The model thus con- 
tains nine free parameters to be determined from fits to the data. 
The best-fit parameters can then be related to physical quantities 
in the system (see, e.g., Sect. 2.2 for the presupernova mass- 
loss rate). Other attenuation processes were taken into account 
for some SNe, e.g. FFA by an ionized medium on the light of 
sight far from the SN progenitor, but there are not important for 
SN 1993J (Weiler et al. [2007l . 

We show here that an estimate of the mean magnetic field 
strength in the radio emission region, (B), can be obtained from 
a comparison of this parametric model with the SSA model of 
Chevalier ( 1998 ). In the latest formalism, the synchrotron emis- 
sion is produced by relativistic electrons having a power-law en- 
ergy distribution: N e (E) = NqE~ 7 , where the spectral index y 
is independent of time. The emitting spherical shell is approx- 
imated by a planar region in the plane of the sky with an area 
nR 2 , where R is the radius. In the absence of FFA by an external 
medium, the observable flux density can be written as 



F v = S v [l - exp(-rssA)] , 
with 

D l c 6 (y)c' 9 (y) \2ci 



5/2 



(A.8) 



(A.9) 



where D is the source distance, the constant c\ = 6.26 x 10 18 in 
cgs units, 05(7), C(,(y), and cg(y) are tabulated as a function of y 
by Pacholczyk ( 1970 ) and 



c' 9 (y) 



V^ F[(y + 6)/4] 
2 r[(y + 8)/4] 



(A. 10) 



T being the gamma function. The factors cg(y) and c' 9 (y) arise 
from the averaging of the synchrotron emission and absorption 



Table A.l. Parameter c(y) for use in Eq. dA.12| ) 
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1.5 


2 


2.5 


3 


3.5 


4 


4.5 


c(y) 


0.724 


0.492 


0.367 


0.289 


0.236 


0.197 


0.169 



coefficients, respectively, over the particle pitch angle < 
1 994 Q. We note that 



C5 ^ 7 ] C f y ] = 3.49 x 10 17 c(y) (cgs units), 
Q,(y)c 9 (y) 



(Longair 



(A.ll) 



where the expression for c(y) can be found in Fransson & 
Bjornsson d 19981 Eq. 29). Calculated values of c(y) are given 
in Table I A.ll for seven values of y. By equating the flux density 
from Eq. ( IA.8b with the one obtained from the model of Weiler 
et al. (Eq.lATT) for AjL f j 
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where we have identified the radius R with the forward shock 
radius R s = R (t/l d ay) m . 

Equation (|A. 121 > can in principle provide an estimate of the 
mean magnetic field in the synchrotron-emitting region as a 
function of the fitted parameters K\, K$, B, and 5" . But in the 
case of combined SSA and FFA, the result strongly depends 
on the FFA model, which is uncertain (Sect. 2.2). Thus, with 
the best-fit parameters B = -0.73 and 8" = -2.05 obtained 
for SN 1993 J by Weiler et al. (120071 ), one gets the unreason- 
able time dependence (B) oc f - 68 (for m = 0.83). We argue in 
Sect. 4.3 that magnetic field amplification by cosmic-ray stream- 
ing should lead to (B) oc t with b ~ 1, as long as the shock is 
not strongly modified by the back pressure from the energetic 
ions. Thus, we expect the relation 



B-8" 



4m + 1 



(A. 13) 



to hold in most of the cases. Using this constraint in the fitting 
procedure could allow to better determine the physical parame- 
ters of radio SNe, in particular the mass loss rate of the progeni- 
tor star. 



Appendix B: Radiative transfer calculations 

The synchrotron emission coefficient (in erg crrT 3 s~' sr~' Hz -1 ) 
averaged over the pitch angle is given as a function of radius by 
(see, e.g., Pacholcz vk l 1 9701 ) 

e v (R) = — I s in <pd</> I p v (E,4>,R)N e (E,R)dE (B.l) 
8?r Jo Jo 

and the absorption coefficient (crrT 1 ) by 

„2 
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7 These quantities were not taken into account by Chevalier J19981 , 
who made the approximation of replacing (B) sin <f> by (B). 
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where p v (E,<p,R) is the total emitted power per frequency per 
electron given in cgs units by the well-known formula 

V3e 3 B(fi)sin0 v f°° 
p v (E,cf>,R)= y K^drj. (B.3) 

Here, -e is the electron charge, K5/3 the modified Bessel func- 
tion, and 

v c = -2Z—B(K)E 2 sia</> (B.4) 
Ajim e c^ 

the critical frequency. In the computer code I have developed, 
the emission coefficient is calculated by numerical integration of 
Eq. ( IB. 11 ) using the tabulation of Pacholczyk (1 19701 1 for the inte- 
gral of the modified Bessel function. The absorption coefficient 
is obtained by approximating the electron energy distribution by 
a power law N e (E c , R) = No(R)E^' / near the characteristic energy 

/ Anmlc'v \ 1/2 

E C = L . J ■ (B.5) 
\3eB(R) sm cf>) 

The absorption coefficient can then be simplified to (Pacholczyk 
[19701 Longair[l99l 

/ v \-<r+4)/2 

k v (R) = c 6 (y)c'g(y)No(R)B(R) (y+2),2 i — \ . (B.6) 

This approximation is used to speed up the computation time. Its 
accuracy has been checked by solving numerically Eq. ( IB. 21 ). 

The brightness profile, I v (p), where p is the distance to the 
center of the disk resulting from the projection of the SNR 
onto the plane of the sky, is calculated from the source function 
S v = £ v /k v by solving the standard transfer equation (Rybicki & 
Lightman fT979l ): 



hip) 



= f e-^Sy(T'y)dT'y, (B.7) 

Jo 
the 
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where the optical depths to SSA are defined as 

-{ 

t v = I K v (R)dz, (B.8) 

f K v (R)dz , (B.9) 

J -(xH(p-p abs ) 



with I = y/R 2 s - p 2 and z = y/R 2 - p 2 . The term H(p - p abs ) 
accounts for attenuation of the emission from the side of the shell 
moving away from us by FFA in the interior of the SNR (e.g. 
Chevalier 1982b); here H denotes the Heaviside step-function 
(H{x) = for x < 0, H(x) = 1 for x > 0). For simplicity, this 
absorption is modeled by a completely opaque disk of radius 
p a bs in the center of the remnant. Other absorption models were 
studied by Bietenholz et al. ( 12003ft . We note that FFA of the 
radio waves in the presupernova wind is not taken into account 
in Eq. jB.lj . In the limit of negligible SSA (i.e. optically-thin 
emission) Eq. ( IB. 71 ) reduces to 



I v (p)= f € v (R)dz. (B.10) 

J-exH(p-p^) 

Finally, the total flux density is given by 



2n C R < 
D 2 i 



pl v (p)dp. (B.ll) 

This quantity is multiplied by the attenuation factor A^^ ps 
(Eq. IA.3b for the modeling of the radio emission at early epochs 
when external FFA is important. 



